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Abstract 

Neutrinos are unique cosmic messengers. Present attempts are directed to extend the win- 
dow of cosmic neutrino observation from low energies (Sun, supernovae) to much higher 
energies. The aim is to study the most violent processes in the Universe which acceler- 
ate charged particles to highest energies, far beyond the reach of laboratory experiments 
on Earth. These processes must be accompanied by the emission of neutrinos. Neutrinos 
are electrically neutral and interact only weakly with ordinary matter; they thus propagate 
through the Universe without absorption or deflection, pointing back to their origin. Their 
feeble interaction, however, makes them extremely difficult to detect. The years 2008-2010 
have witnessed remarkable steps in developing high energy neutrino telescopes. In 2010, 
the cubic-kilometre neutrino telescope IceCube at the South Pole has been completed. In 
the Mediterranean Sea the first-generation neutrino telescope ANTARES takes data since 
2008, and efforts are directed towards KM3NeT, a telescope on the scale of several cubic 
kilometres. The next years will be key years for opening the neutrino window to the high 
energy Universe. With an instrumented volume of a cubic kilometre, IceCube is entering a 
region with realistic discovery potential. Discoveries or non-discoveries of IceCube will have 
a strong impact on the future of the field and possibly mark a "moment of truth". In this 
review, we discuss the scientific case for neutrino telescopes, describe the detection princi- 
ple and its implementation in first- and second-generation installations and finally collect 
the existing physics results and the expectations for future detectors. We conclude with 
an outlook to alternative detection methods, in particular for neutrinos of extremely high 
energies. 
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1. Introduction 



High-energy neutrinos, witli energies much larger than 100 MeV, must be emitted as a 
by-product of colhsions of charged cosmic rays with matter; in fact, only neutrinos provide 
incontrovertible evidence for hadronic acceleration. Since they can escape much denser 
celestial environments than light, they can be tracers of processes which stay hidden to 
traditional astronomy. At the same time, however, their extremely low reaction probability 
makes their detection extraordinarily difficult. 

First ideas to search for cosmic neutrinos other than those from the Sun date back to 
the late fifties. In 1960, K. Greisen proposed a 3000 ton underground Cherenkov detector 
to record neutrinos emitted by the Crab nebula [IJ. He was fianked by F. Reines pj, who 
realised, however, that "the cosmic neutrino fiux cannot be usefully predicted" - i.e. that 
a mass of 3000 tons may be far too small. In the same year, M. Markov made his ground- 
breaking proposal "to install detectors deep in a lake or in the sea to determine the direction 
of charged particles with the help of Cherenkov radiation" p3]. In the decades since then it 
was realised that high-energy neutrino astronomy requires detectors of a cubic kilometre or 
larger that can indeed only be implemented in open media. Actually, the first project of 
that size, IceCube at the South Pole, has just been completed. Two others, KM3NeT in the 
Mediterranean Sea and GVD in Lake Baikal, are in their preparatory phases. No doubt, we 
are entering an exciting era of opportunity. 

Already now, neutrino astronomy is reality in the low- energy sector, where the detection 
of neutrinos from the Sun and the supernova SN 1987A was honoured by the 2002 Nobel Prize 
for physics. Figure [T] shows a compilation of the spectra of dominant natural and artificial 
neutrino fiuxes. No practicable idea exists on how to measure the neutrinos of the 1.9 K 
neutrino counterpart to the cosmic microwave background. At higher energies, neutrinos 
from the Sun, from SN 1987A, from reactors and from the interior of the Earth have already 
been detected, as have so-called "atmospheric neutrinos" created in cosmic ray interactions 
in the Earth's atmosphere. Still awaiting detection are high-energy cosmic neutrinos from 
extraterrestrial sources such as active galactic nuclei (AGN) or from interactions of ultra- 
energetic protons with the cosmic microwave background [1]. These cosmic neutrinos will 
hopefully be detected by neutrino telescopes in the next decade, even though predictions for 
their fiuxes are uncertain by orders of magnitude in many cases. 

The development of high-energy neutrino astronomy is refiected in a series of previous 
reviews spanning the period 1995 to 2009 [51-ITO]. The neutrino telescopes discussed in this 
review focus on energies beyond a few GeV. First searches for such neutrinos were made 
in the 1960s in the Kolar Gold Field mine in India and in the East Rand mine in South 
Africa (for a review see [7|). In the 1980s, the spectrum of atmospheric muon neutrinos was 
measured with a detector in the Frejus tunnel between France and Italy, and a first limit on 
the diffuse flux of extra-terrestrial TeV neutrinos was set [11]. Over the following decades, 
the evolution of underground neutrino detectors culminated in two experiments with an 
area of about 1000 m^ each (see Sect. 3.5 for a discussion of effective areas): MACRO in 
the Gran Sasso Underground Laboratory in Italy and Super-Kamiokande in the Japanese 
Kamioka mine. MACRO collected more than thousand atmospheric neutrinos over six years 
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Figure 1. Measured and expected 



Neutrino energy fluxes of natural and reactor neutrinos. 

of data taking. Super-Kamiokande, with an even larger data sample, is still in operation. 
The atmospheric neutrino results from these detectors have demonstrated that neutrinos 
oscillate between their flavour states and u^, additionally to the oscillations observed 
for solar neutrinos [7]. 

The first-generation detectors in water and ice have beaten the largest underground 
detectors by a factor of about 30 with respect to their sensitivity to high-energy neutrinos. 
The second-stage detectors on the cubic-kilometre scale will yield another factor of 30. 
Compared to detectors underground we therefore enter a "factor-1000 era". Arguably, this 
factor is not a guarantee for discoveries. On the other hand it rarely happened in astronomy 
that improvements of more than an order of magnitude (in sensitivity or in angular or time 
resolution) came along without discovering new, unexpected phenomena |12j. "Nothing is 
guaranteed, but history is on our side" [13]: In some years we will know whether we indeed 
have entered an era of discovery or not. 

This review is organised as follows: Section |2] summarises the scientific motivation. Apart 
from the main topic, neutrino astrophysics, it includes the indirect search for dark matter, 
the study of standard and non-standard neutrino oscillations, the search for exotic particles 
like magnetic monopoles, super-symmetric Q-balls or nuclearites and - last but not least - 
the investigation of environmental effects, be it in deep natural water or Antarctic ice. The 
basics of the detection methods are summarised in Sect. |3} In Sect. |4]the first-generation 
neutrino telescopes are described, in Sect. [5] the second-generation projects on the cubic- 
kilometre scale. A selection of results obtained with NT200 in Lake Baikal, ANTARES in 
the Mediterranean Sea as well as AMANDA and IceCube at the South Pole is presented 
in the following Sect. [6j For the highest energies beyond 100 PeV, even cubic-kilometre 
detectors are far too small to detect the feeble neutrino fluxes expected. This is the realm 
of new technologies which aim, with a correspondingly high detection threshold, to monitor 
volumes of 100 cubic kilometres and beyond. These methods are described in Sect. [7j The 
last section finally gives a summary and tries an outlook to forthcoming developments. 



2 



2. Scientific Background and Motivation 



The primary motivation to build kilometre-scale neutrino detectors is driven by the ob- 
servation of charged cosmic rays. Since long, neutrinos have been supposed to be a key 
messenger to identify the sources of cosmic ray acceleration and to provide a deeper under- 



standing of the associated astrophysical objects. In Sect. 2.1, we recall basic information on 
cosmic rays and gamma rays which motivate (and constrain) the search for energetic neutri- 
nos. We sketch the astrophysical source candidates for neutrinos from hadronic acceleration 
processes and relate the predicted fluxes to event rates expected in neutrino telescopes on 
the cubic-kilometre scale. 

The identification of the sources of cosmic rays and, more generally, the opening of a 
new observational window to the Universe, is arguably the most important, but by far not 
the only purpose of large neutrino telescopes. Neutrino detectors are multi-purpose devices 



addressing also questions of particle physics and environmental science. Section 2.2 gives 



an overview on the particle physics issues studied with the help of neutrino telescopes, and 



Sect. 2.3 sketches the environmental spin-offs. 



2.1. Neutrinos from cosmic accelerators 

High-energy charged cosmic rays have been known to exist for almost a century, but their 
origin is still a mystery. Since charged particles are deflected in the inter- and extragalactic 
magnetic fields, their arrival direction at Earth does not reveal their sources, except at the 
very highest energies. In contrast, neutrinos, produced at the acceleration sites or during 
cosmic ray propagation, propagate on straight trajectories and point back to their origin. 

Several astrophysical object classes have been proposed as potential particle accelerators. 
In spite of a vast amount of observational data in all electromagnetic wavelength regimes - 
from radio to gamma - it is, however, still unclear, whether the non-thermal processes in 
these objects are of electronic or hadronic nature. It appears likely that a final answer to 



this question will require the observation of neutrinos (cf. Sect. 2.1.2). 



2.1.1. Cosmic rays 

Cosmic rays are charged particles - basically hadrons (protons, light and heavy nuclei) 
with only a tiny admixture of electrons - that impinge on Earth from outer space, with an 
energy spectrum extending to energies in excess of 10^° eV (100 EeV), see Fig. [2} 

The energy spectrum follows a broken power law E'"". The region where the decrease 
with a power index a = 2.7 steepens to a power index of a = 3.1 at about 10^^ eV (1 PeV) 
is called the "knee". It is assumed that cosmic rays up to (and even beyond) the knee are 
of Galactic origin. The main accelerator candidates are supernova remnants (SNR): The 
shock fronts powered by supernova explosions propagate into the interstellar medium, and 
by repeated scattering processes across the shock front particles can gain energy ("first-order 
Fermi acceleration" [T514T9] ). 

The spectrum resulting from first-order Fermi acceleration can be shown to follow roughly 
an spectrum, with a maximum energy 

Emax (X Z ■ V B ■ L , (1) 



3 



air shower measurements 



direct measurements 



\ 



■ direct measurements (protons) 
• direct measurements (all particles) 
A air siiower data (all particles) 



I 











10 



10 



10 10 



10 



10 



10 



10 10 



10 



10 



Figure 2. Energy spectrum of 
cosmic rays, illustrating its main 
features (knee, ankle, and 
high-energy cut-off at a few 
10^^ GeV). Figure taken from 



Energy En [GeV] [JJ]. 



where Z is the charge of the particle, v the velocity of the shock wave, B the magnetic field 
strength in the acceleration region and L the size of this region. 

With eq. g, the maximum energy to which a SNR can boost particles turns out to be 
somewhat below the knee, at about 10^^ eV. However, by interactions of cosmic rays with the 
magnetic fields in the acceleration region, these fields can be amplified, resulting in energies 
of up to 10^^ eV (see [20] and references therein). The mechanism of converting kinetic 
energy of the SNR into energy of accelerated particles is effective over the first 10^-10^ years 
of the SNR. 

Circumstantial evidence for SNRs as the main sources of cosmic rays can be also derived 
from the Galactic energy budget. The energy density of Galactic cosmic rays is about 
leV/cm^ [21j, corresponding to pcR = 10~^^ erg/cm^. The power required to sustain this 
density is L = pcR ■ ^/^cont? with tcont ~ 10^ years [14j being the average containment 
time of cosmic rays in the Galaxy and V ~ 10®^ cm^ the volume of the visible part of the 
Galaxy. With L ?a 10^^ erg/s, the power to maintain the cosmic ray density turns out to be 
about 10% of the power generated by supernovae, which release on average 10^^ erg every 
20-50 years, a relation noticed already in 1934 by Baade and Zwicky |22| . 

Besides supernova remnants, there are further candidates for Galactic particle accelera- 
tion, most notably pulsars with their extremely high magnetic fields at the poles, and binary 
systems with a neutron star or a black hole as one of the partners. The latter often form 
relativistic radio jets and are then dubbed microquasars (for a review, see [23j). 

Somewhere between 10^^ and 10^^'^ eV, at maximum, known Galactic source candidates 
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are running out of power and extragalactic sources start dominating the spectrum. Around 
j^gis.s gY j^YiQ spectrum flattens again to an E^'^''^ shape, a feature named the "ankle". Almost 
all cosmic rays with energies above the ankle are assumed to be of extragalactic origin. Even 
if Galactic sources would accelerate particles beyond the ankle energy, these particles would 
right away escape the Galaxy since their gyroradius exceeds the size of the Galaxy. 

The cut-off of the cosmic ray spectrum at highest energies, confirmed by recent Auger 
measurements is likely due to interactions with the cosmic microwave background 

{p + 7cMB — ^ and was first predicted by Greisen, Zatsepin and Kuzmin [25l |26] : the 
"GZK cut-off". 

The main candidates for particle acceleration to energies beyond 10^^ eV are Active 
Galaxies, Gamma Ray Bursts and starburst galaxies. With eq. ([T]), the maximum energy to 
which these objects could accelerate protons is between 10^*^ and 10^^ eV. Similarly as for 
supernova remnants and the energy budget of the Galaxy, there is a suggestive coincidence 
between the power released by these extragalactic objects and the power required to sustain 
an E'"^ flux up to GZK energies. The latter is 10"'^ eV/cm'^, equivalent to 3 x 10~^^ erg/cm^ 
|27| . Normalised to the cosmic abundance of AGN and GRB, this density translates to a 
required power of about 2 x lO^^erg/s (Active Galaxies) and 3 x lO^^erg/s (Gamma Ray 
Bursts) and this indeed is of a similar order as the corresponding electromagnetic energies 
released by these objects. 

For a recent review of cosmic rays we refer to [14j). 

2.1.2. Production of neutrinos 

It has been recognised for half a century [H |2] that protons from cosmic accelerators 
would also generate neutrinos, via charged pion production in collisions with the ambient 
matter or radiation fields, in reactions such as: 



p -I- nucleus -)■ tt + X (vr = 7r=^, vr") (2) 
P + 7 ^ A+^<^ ^ / , (3) 



with the subsequent decays tt — )■ 77, tt — )■ ^ and — e^V^Vf,, /i~ — e'v^Vf.. The 
resulting neutrino flavour ratio is approximately z/g : z/^ : i^,- = 1 : 2 : at the sources; 
neutrino oscillation turns this into a ratio of z/g : z/^ : z^r = 1 : 1 : 1 upon arrival at Earth 



(see Sect. 2.2.3). 



The kinematic threshold for process ^ is determined by the photon energies in the 
radiation field. For ambient photons in the UV region, as characteristic for many stars and 
accreting objects, it is in the range of several PeV. For the cosmic microwave background, it 
is at about 10^^'^ eV. If the photon spectrum has a broad spectrum, such as the power-law 
for photons generated by synchrotron radiation, the threshold is "smeared" to much lower 
energies. This is particularly important for neutrino astronomy in the TeV range. 

TeV gamma rays can be produced via the decay of neutral pions, but also by inverse 
Compton scattering: 

^ ~l~ 7low energy ^ *^ ~l~ 7high energy ■ (4) 
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Actually, most of the measured spectra from TeV gamma ray sources are compatible with 
models based on inverse Compton scattering, and in many cases with the so-called syn- 
chrotron-self Compton (SSC) model where the photon gas is provided by synchrotron radia- 
tion from accelerated electrons (see Fig.js]). Needless to emphasise that pure SSC models are 
based on electron, not hadron acceleration and do not directly explain the origin of cosmic 
rays. Pure electron acceleration models are called "leptonic models". In most realistic cases, 
both electrons and hadrons will be accelerated. 

Figure |3] sketches the processes happening in such a combined model. The synchrotron 
radiation from electrons serves as target for Inverse Compton scattering as well as for proton 
collisions. Electrons are cooled by synchrotron emission and may boost synchrotron photons 
to the 10-100 TeV range but certainly not to PeV energies. The observation of PeV gamma 
rays would therefore be a clear proof of hadron acceleration. Unfortunately, the range of 
PeV photons does not exceed the size of our Galaxy, since they are absorbed by the process 
7PeV + 7CMB e+e". 



high en. 




Figure 3. Gamma ray and neutrino production in a jet 
emitted from an Active Galactic Nucleus, with both hadrons 
and electrons being accelerated along the jet. 



If charged particles are confined by large magnetic fields, only neutral particles can escape 
the acceleration region. Apart from gamma rays and neutrinos, this can be also neutrons 
(see top right in Fig. |3|, provided the source is sufficiently dilute. Neutrons decaying outside 
the source would then yield those cosmic ray protons which are detected on Earth and which 



can be used to constrain the flux of neutrinos (see Sect. 2.1.4). 



2.1.3. Galactic sources 

Until the mid-1990s, only one supernova remnant, the Crab nebula, had been detected 
in TeV gamma rays. Therefore, predictions for neutrinos from these sources were not yet on 
a firm ground at that time. The Crab belongs to a special class of SNR, the pulsar wind 
nebulae, where a central pulsar emits material into the nebula. The most likely sources of 
Galactic cosmic rays are shell-type supernova where this effect is absent or not essential. 
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In the mean time the situation has changed dramatically. Imaging Air Cherenkov Tele- 
scopes like H.E.S.S., MAGIC and VERITAS have detected more than hundred sources of 
TeV gamma rays, amongst them about 30 SNRs (for recent reviews see [281 129]). From the 
observed fluxes of gamma rays, estimates or upper bounds on the flux of neutrinos can be 
derived, assuming that most of the observed gamma rays stem from decays of tt^s generated 
according to eqs. ^ and (|3]). Keeping in mind that high-energy gamma rays may well 
emerge from inverse Compton scattering, eq. (|4]), their observation is not a proof that the 
source accelerates hadrons rather than only electrons. A certain test can be provided by 
detailed information on the MeV-GeV part of the gamma spectrum and by information on 
its high-energy cut-off. For Galactic sources, the morphology of gamma ray emission can 
be studied and provides additional information. In this context, we note the observation 
of the supernova remnant RXJ1713. 7-3946 with the H.E.S.S. telescope [30J. The image 
of this source (see Fig. |4| left) shows an increase of the gamma flux from the direction of 
known molecular clouds. The effect can be attributed to protons accelerated in the SNR 
and then interacting with the clouds. The spectrum of gamma rays above some 100 GeV 
is well compatible with expectations for the decay of 7r°s from proton interactions. Recent 
measurements with the Fermi-LAT instrument |31] . however, indicate that the spectrum 
at lower energies may be better described by leptonic models. This is in accordance with 
X-ray line emission around 1 keV which should be produced along with pion decays but has 
not been observed with the Suzaku satellite detector [32j. Unambiguous Galactic sites of 
hadronic acceleration have thus yet to be identified. Extended gamma ray sources which 
trace the density of molecular clouds, such as the SNR Vela (RX J0852. 0-4622) and a 
large region near the Galactic Centre |34| . remain good candidates. 




Figure 4. Left: Gamma ray image of the supernova remnant RX J1713. 7-3946, recorded with 
the H.E.S.S. telescope; superimposed are contours of X-ray surface brightness recorded by the 
ASCA satehite (figure taken from |35]). Right: Measured gamma flux from RXJ1713. 7-3946 and 
the related neutrino flux estimated under the assumption that the gamma emission is of purely 
hadronic origin (figure taken from [36])- 
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In |36| . the expected neutrino flux from tlie SNR RX J1713. 7-3946 is calculated (cf. 
Figj4| right). Based on the assumption that all observed gamma rays stem from vr" decays, 
the authors calculate the neutrino flux expected from vr^ decays. For five years of data 
taking with a cubic-kilometre detector and a threshold at 1 TeV, the number of events is 
calculated to be between 7 and 14, over a background of 21 atmospheric neutrino events in 
a 2.1° search cone (resulting from the 1.3° diameter of the source, multiplied with a factor 
of 1.6 to achieve optimal sensitivity). A threshold of 5 TeV results in 2.6-6.7 signal events 
and 8.2 atmospheric neutrino events. Cutting at higher energies will eliminate not only 
the atmospheric background but also the signal, cutting at lower energies will significantly 
worsen the signal-to-background ratio. These estimates suggest that the positive effect of 
a detector threshold much below 1 TeV will be small, at least for steady sources. It is also 
clear that neutrinos from this particular source will be hardly observable if the gamma flux 
of hadronic origin is sub-dominant; in this case, however, the quest of identifying the sites 
of hadronic acceleration starts afresh, and neutrino observations will be more crucial than 
ever in solving this puzzle. 

Several candidates for sources with hadron acceleration beyond the knee ("Pevatrons") 
have been identified in the Cygnus region by the Milagro collaboration |H7 t [55]. The gamma 
ray spectrum of the strongest of these sources, MGRO J1908+06, is consistent with an E~'^ 
behaviour between 500 GeV and 40 TeV. It does not show evidence for a cut-off, in accor- 
dance with the expectations for Pevatrons, which should emit gamma rays with energies 
up to several hundred TeV. In jES], the associated neutrino fluxes from the six identified 
Milagro source candidates have been calculated and the event rates in the IceCube neutrino 
telescope estimated. In a simulated neutrino sky-map for 5 years of data taking, two of the 
sources are discernible "by eye" when applying a lower energy threshold of 40 TeV. Stacking 
all six sources, an excess is found with a Poisson probability for being a background fluctu- 
ation of smaller than 10"'^ for lower energy thresholds anywhere between 10 and 100 TeV. 
The simulation assumes a gamma-ray cut-off at 300 TeV. 

Similarly close to the sensitivity of cubic-kilometre telescopes are the expected neutrino 
fluxes from microquasars, as for instance LSI +61 303 and LS5039, where the latter could 
provide a handful of neutrino events per year |40J. 

Practically all publications of the last years (see references above and e.g. |36| HTI - 145) ) 
come to the unanimous conclusion that cubic-kilometre detectors will just "scrape" the detec- 
tion region. The present estimates suggest that the sensitivity of a cubic-kilometre telescope 
is "tantalisingly (and frustratingly) close" to the expectations for the brightest observed 
Galactic TeV gamma sources. 

2.1.4- Extragalactic sources 

There is much less observational guidance to predict neutrino fluxes from extragalactic 
than for Galactic sources. The predictions for individual sources have order-of-magnitude 
character. The best-motivated extragalactic candidates for high energy neutrino emission 
are Active Galactic Nuclei (AGN), Gamma Ray Bursts (GRB) and galaxy clusters. 

• Active Galactic Nuclei host a super-massive black hole with 10^-10^ solar masses in 
their extremely bright centre. The black hole accretes matter and thus transforms huge 
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amounts of gravitational energy into radiation, typically of the order of 10^^ erg/s. This 
energy can also be converted to kinetic energy of accelerated particles, see also Fig. |3j 
Blazars, a particular class of AGN where the jet is aligned closely to the line of sight, 
turned out to be strong gamma emitters, with 119 sources at GeV energies listed in 
the 2009 Fermi bright gamma source list [46j and 18 at TeV energies being reported in 
|47| . Gamma ray emission from blazars is often highly variable, with the most extreme 
variation observed by H.E.S.S. for the blazar PKS 2155-304: An increase by two orders 
of magnitude within one hour [48j . Naturally, the observation of neutrino events from 
such a source and within such a short time window would be rather significant. Another 
interesting outburst has been observed from the blazar lES 1959+650 ^49j. This was 
an "orphan fiare", where the TeV emission was not accompanied by X-ray emission, 
as it typically would be for SSC models. A hadronic model does not require such a 
correlation between TeV and X-ray (synchrotron) emission; therefore orphan fiares are 
interesting environments to search for neutrino emission. 

• Gamma Ray Bursts are the most cataclysmic phenomena in the Universe, releasing 
huge amounts of energy in gamma rays within milliseconds to minutes. The favoured 
explanation for the longer bursts is the collapse of a massive star into a black hole. The 
so-called fireball model of GRB assumes that a central engine ejects large amounts of 
mass within a short time interval which form successive plasma shells and have typical 
Lorentz factors of F = 100-1000. When the outer shells slow down they are hit by 
the inner shells and internal shock fronts are piling up. Along these fronts, electrons 
and protons are accelerated. Electrons are cooled by synchrotron radiation, protons 
can be accelerated up to energies as high as 10^^ eV |i50j. When the shells run into the 
interstellar medium, external shocks are built up, leading to afterglow emission visible 
in X-ray, optical and radio wavelengths. Neutrino emission has been calculated for 
three phases: the precursor phase when the jet is still forming and no electromagnetic 
radiation is escaping [51j; the "prompt" phase coinciding with the burst in gamma rays 
(see e.g. |,50j); and the afterglow phase [52j. We will return to these predictions when 
presenting experimental bounds in Sect |6| 

• Starhurst galaxies are galaxies undergoing an episode of large-scale star formation, 
where the central regions eject a galactic-scale wind driven by the collective effect of 
supernova explosions and winds from massive stars. Recently, the starburst galaxies 
NGC253 (southern hemisphere) and M82 (northern hemisphere) have been detected 
by the H.E.S.S. and VERITAS telescopes IM] . The gamma ray flux at several 
hundred GeV suggests cosmic ray densities of two to three orders of magnitude above 
that in our own Galaxy. Following [55j, the cumulative neutrino flux of all starburst 
galaxies is detectable with cubic-kilometre detectors. 

Predictions for the integrated flux from all extragalactic sources are based on the observed 
gamma and X-rays fluxes or those of charged cosmic rays above 10^^ eV. A selection of 
predictions and bounds is shown in Fig. [s] Early normalisations as e.g. in [56j assumed 
neutrinos would be produced in the cores of AGN, accompanied by X-ray emission, and 
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that most of the observed X-ray background was non-thermal radiation from a superposition 
of the fluxes from unresolved AGN. This model as well as others violated upper bounds 
derived from the observed cosmic ray spectrum. Subsequent observations of these AGN 
showed that most of the X-ray emission is thermal and therefore cannot be directly related 
to the production of relativistic particles. Relying on MeV gamma ray rather than X-ray 
observations, the authors scaled down the original prediction by an order of magnitude 
cf . curve 2 in Fig. [s]) . 

Figure 5. Spectra of diffuse neutrino fluxes, 
multiplied with E"^. Points: AMANDA 
measurements of the atmospheric neutrino flux 
[58], compared to the prediction from |59) 
(dotted line). Horizontal band: E~'^ upper 
bound derived from charged cosmic rays at 
E > 10^^ eV [60], with the width reflecting the 
uncertainty of cosmic evolution parameters. 
Numbered curves: predictions for the diffuse 
flux of neutrinos from all AGN (1,2) [571 E], of 
cosmogenic neutrinos (3) [3] and of neutrinos 
from GRBs, assuming that they are the sources 
of highest energy charged cosmic rays (4) |50) . 
Figure courtesy of J. Becker. 

Two upper bounds on the diffuse neutrino flux, both derived from charged cosmic ray 
fluxes, are frequently used as benchmarks. The first ("Waxman-Bahcall bound" [60j) is 
normalised to the cosmic ray flux at about 10^^ eV. Assuming a generic spectrum for all 
extragalactic sources, the authors obtain a limit of E"^^^ = 1-5 x 10^^ GeV cm^^ s^^ sr~^, 
where is the neutrino flux differential in energy, time, area and solid angle and the 
uncertainty is given by different cosmic evolution models (coloured band in Fig. [s]). This 
estimate assumes that the sources are sufficiently dilute, so that neutrons can escape, decay 
and provide the observed cosmic rays to which the estimate is normalised. If the sources are 
opaque even for neutrons, the only remaining estimator is electromagnetic radiation. The 
decay of tt^s (co-produced with charged pions producing neutrinos) yields gamma rays which 
develop electromagnetic cascades. The energy of gammas escaping the source is mostly in the 
range 1 MeV to 100 GeV. The diffuse gamma-ray background above 30 MeV was measured 
by the EGRET satehite [62j to be E^*^ = 1.37 x 10"^ GeV cm'^ s"^ sr'^ and sets a bound of 
similar size to the neutrino flux (not shown in the figure), which hardly can be circumvented 
by more sophisticated assumptions on the character of the sources. Recent Fermi data |63] 
suggest [64] that this "gamma bound" (or "cascade bound", originally derived in [65|) is close 
to the Waxman-Bahcall bound, and actually has already been superseded by experimental 
upper limits from the running underwater/ice telescopes (see below). In a cubic-kilometre 
detector, the Waxman-Bahcall flux would lead to 100-500 events per year. 

Contrary to Waxman and Bahcall, Mannheim, Protheroe and Rachen (MPR) |61) as- 
sume that a significant part of the observed cosmic ray spectrum between 10^^ eV and 
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10^^ eV is due to extragalactic rather than Galactic sources. Interpreting this spectrum 
as a superposition of spectra from many extragalactic source classes, each with a differ- 
ent cut-off, the neutrino bound considerably weakens towards lower energies and is about 
E'^ ■ dN/dE ^ 5 X 10""^ GeVcm"^ s~^ sr~\ at a few 10^"^ eV. Meanwhile, also this bound is 
only of historical interest, due to the new results from the Fermi satellite and from neutrino 
telescopes. 

A non-astrophysical contribution to the diffuse neutrino background is assumed in "top- 
down scenarios", where cosmic rays of highest energy are due to cascade decays of long-lived 
super-heavy particles at unification-scale energies (10^^-lO^^eV). Such processes inevitably 
produce high-energy neutrinos and photons. An estimate of the resulting neutrino event 
rate in IceCube of up to 40 events per year |66| has meanwhile been superseded by recent, 
tight limits from Auger on the photon contribution to energetic cosmic rays |67] . 

The most exciting discovery of neutrino astronomy would be the detection of point 
sources rather than just a high-energy excess in diffuse fluxes. However, experimental limits 
for diffuse fluxes are setting bounds for expected point-source fluxes. The argument is as 
follows: Contributions to the diffuse flux will come from all the observable universe, up to a 
distance c/Hq, whereas point sources will be visible at a level of several events per source only 
up to a limited distance of a few hundred Mpc, assuming reasonable maximum luminosities 
per source. For a homogeneous distribution of extragalactic sources, one therefore can derive 
a limit on the number of observable point sources. In jH] the following assumptions are 
made: a homogeneous source density in a Euclidean universe; a "typical" (and similar) 
source luminosity Lgource for all sources; an E'"^ spectrum of the neutrino fluxes. Given 
an experimental limit -ft'diffuse on the diffuse neutrino flux and a sensitivity Cpoint to point 
sources, the expected number of resolvable extragalactic point sources, A^^ scales as 

7> r -^diffuse source / ^ \ 

^ • (5) 

point 

With the present diffuse flux limit from IceCube and its point source sensitivity, one obtains 
Ns ~ 0.2-2 (in analogy to an estimate made in |68) for the AMANDA diffuse limit). This 
means that, under the given assumptions, a cubic-kilometre detector would have a remaining 
chance to detect extragalactic point sources. Note, however, that a few individual, very close 
sources could circumvent the homogeneity assumption and be well observable. Also, point 
sources with cut-offs below a few hundred TeV would not be covered by the argument above 
since, in order to obtain the best sensitivity for diffuse fluxes, a lower energy cut at about 
100 TeV has to be used [69]. In particular this implies that the argument of [?TJ does not 
apply to Galactic sources since they are not homogeneously distributed and typically are 
expected to have cut-off energies below 100 TeV. 

2.2. Particle physics and exotic phenomena 

Beyond charting the high energy universe, neutrino telescopes may reveal first signatures 
of new physics beyond the standard model. Opportunities include the indirect search for 
dark matter particle candidates, the search for super-heavy exotic particles like supersym- 
metric Q-balls or magnetic monopoles, or the search for deviations from the established 
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neutrino oscillations which may result from violation of Lorentz invariance. We sketch these 
opportunities in the following. 

2.2.1. Indirect search for dark matter 

Favoured dark matter candidates are Weakly Interacting Massive Particles (WIMPs). 
They are preferentially discussed in the framework of the minimal supersymmetric standard 
model (MSSM), see [70j and references therein. Specifically, the lightest supersymmetric 
particle is considered a good WIMP candidate, with predicted masses in the range from a 
few GeV to a few tens of TeV jTT]. A more constrained class of theories is based on minimal 
supergravity (mSUGRA); see [72] for a comparison of MSSM and mSUGRA. WIMPs can 
be trapped in the gravitational potential of celestial bodies such as the Earth or the Sun 
through elastic scattering and would accumulate in the cores of these objects. If WIMPs are 
Majorana particles (as in supersymmetry) , they may subsequently annihilate once sufficient 
densities are reached. Neutrinos are produced via decay of the annihilation products |73] . 
The mean neutrino energy depends on the WIMP mass and the annihilation final state. For 
instance, quark-antiquark states result in softer, W^W~ states in harder neutrino spectra. 

The capture rate depends on the mass of the celestial body and on local density, mean 
velocity and scattering cross section of the WIMPs on nuclei. Slower WIMPs are captured 
more efficiently than faster WIMPs. The scattering cross section can be spin-dependent or 
spin-independent; the capture rate thus depends on the spin composition of the nuclei in 
the celestial body (capture in the Sun, e.g., is dominated by scattering on hydrogen nuclei 
with spin 1/2 and is therefore particularly effective for a spin-dependent interaction). The 
scattering cross section also depends on the mass ratio of the WIMP and the scattered 
nucleus. For instance, the Sun with its lighter nuclei captures light WIMPs more efficiently, 
the Earth gains with respect to heavier WIMPs. The annihilation rate scales with the square 
of the accumulated WIMP density (which likely has reached equilibrium for the Sun but 
less likely for the Earth) and with the annihilation cross section (which is related to the 
scattering cross section through the underlying theoretical model). 

WIMP detection via secondary particles from annihilations (in particular gammas or 
neutrinos) is called "indirect detection", in contrast to "direct detection" of recoil nuclei from 
elastic WIMP-nucleus scattering. Results from both methods can be related (see Sect. |6]), 
but one has to keep in mind that direct event rates are highest for fast WIMPs and scale 
with the WIMP density, whereas indirect rates are highest for slow WIMPs and scale with 
the square of the WIMP density. In this respect, the two methods are complementary. 
Whereas direct and indirect searches are on similar footing with respect to spin-independent 
scattering, indirect searches win when it comes to neutralino capture via spin-dependent 
scattering in the Sun |74) . Consequently, limits on the spin-dependent scattering cross 
section from indirect searches are by nearly two orders of magnitude better than those from 
direct searches. For detailed reviews on dark matter searches, we refer to [7^76) . 

An alternative and in many ways orthogonal approach is the search for WIMP candidates 
at accelerators, in particular at the LHC |77| . The forthcoming years will be particularly 
interesting in this respect. One has to keep in mind, however, that a discovery at LHC 
would not yet pin down the dark matter properties since it would remain to be shown that 
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the LHC particles actually (exclusively) form the astrophysical dark matter. In that sense, 
a full picture of the very nature of dark matter (if it were supersymmetric particles) will be 
obtained only by combining results from all three methods |75] . 

2.2.2. Super-heavy particles 

Neutrino telescopes can be used to search for several exotic objects on the mass scale 
much beyond that of WIMPs which may lead to characteristic signal patterns. Here, we 
mention the three most popular ones: magnetic monopoles, supersymmetric Q-balls and 
nuclearites. The predicted mass of magnetic monopoles ranges between 10^ and 10^^ GeV, 
nuclearites could have masses from a few hundred GeV up to the mass scale of neutron 
stars, and Q-balls up to 10^'' GeV. The expected signals for these particles are similar; a 
clear assignment is therefore excluded for a single event, but may be possible if angular and 
light yield distributions can be studied for a sample of several events. 

• Magnetic monopoles: Magnetic monopoles have been introduced by P. Dirac in 1931 
in order to explain the quantisation of electric charge jTS]. After decades of unsuc- 
cessful monopole searches at accelerators and in cosmic rays, the efforts gained new 
momentum when monopoles turned out to be a consequence of most variants of Grand 
Unified Theories [79j. A phase transition in the early Universe at 10~^^s might have 
filled the Universe with a significant amount of monopoles. 

Measurements and estimates of cosmic magnetic fields suggest that they could accel- 
erate magnetic monopoles lighter than 10^^ eV (sub-GUT-scale) to relativistic veloc- 
ities [sot [8T] . The magnetic charge of monopoles obeys the Dirac quantisation rule 
g = n - e/{2a), with n = 1, 2, 3, . . . and a = 1/137 (e being the elementary charge). A 
magnetic monopole with unit charge g = 137e/2 and a velocity above the Cherenkov 
threshold in water or ice (/3 > 0.75) would emit Cherenkov light along its path. The 
light intensity would exceed that of a minimally ionising muon by a factor of 8300, 
this providing a rather unique signature in a neutrino telescope. 

Typical GUT versions predict monopoles with masses of 10^^ GeV and more. These 
monopoles would have typical virial velocities peaking at /3 = 10~^- 10~^. They might 
catalyse baryon decays along their path which would be detected via the decay particles 
[82j . For certain regions of the parameter space spanned by catalysis cross section and 
monopole velocity, the Cherenkov light from the secondary particles would create a 
pattern of a slowly propagating light source on a straight trajectory, well observable 
in a neutrino telescope. 

• Nuclearites: Nuclearites ("strange quark matter" or "strangelets") are hypothetical 
aggregates of u, d and s quarks combined with electrons, to adjust electric neutrality. 
They might be stable for baryon numbers ranging from those of usual nuclei up to 
those of neutron stars [A ^ 10^^) |83| [M]. Nuclearites could have been produced 
in the primordial Universe or in certain astrophysical processes like the collision of 
neutron stars. They would induce a thermal shock wave along their path through the 
detector medium. For virial velocities, the corresponding Planck radiation could reach 
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10^ K. A light source so slow and bright would produce a very specific time pattern 
similar to that of a slow GUT monopole. 

• Supersymmetric Q-balls: Q-Balls are hypothesised coherent states of squarks, sleptons 
and Higgs-fields |85| . These SUSY soliton states could be stable for masses M ^ 
10^^ GeV, the life time for smaller masses is unknown. Stable Q-Balls could provide 
a relevant contribution to Dark Matter. Decays of unstable Q-balls could produce 
WIMPs, with properties not compatible with the conventional thermal origin [8A{ . 
The discovery of Q-Balls would have an enormous impact on the understanding of 
matter-antimatter symmetry emerging in the early Universe |86]488] . 

Neutral Q-ball objects (SENS, supersymmetric Electrically Neutral Solitons) could 
catalyse proton decays along their path, similar to GUT monopoles. Electrically 
charged Q-ball objects (SECS, supersymmetric electrically charged solitons) would 
produce light in a similar way as nuclearites. 

2.2.3. Neutrino oscillations 

We meanwhile know that neutrinos change their fiavours during propagation. There are 
two standard mechanisms that contribute to this phenomenon, a mixing between mass and 
fiavour eigenstates (inducing "vacuum oscillations"), and fiavour-dependent forward scatter- 
ing amplitudes in matter ("MSW effect" jSHUSO])- Non-standard oscillation effects could e.g. 
be associated to a violation of Lorentz invariance. Neutrino telescopes have the potential to 
study the high-energy regime of neutrino oscillations. 

Standard oscillations: 

Neutrino data from solar, atmospheric, reactor and accelerator experiments support the 
concept of neutrino oscillations. The weak flavour eigenstates I'eji^fi, i't are linear combina- 
tions of mass eigenstates z/i,z/2,t'3. For the simplifled case of two flavours z/^, i^r and two 
mass eigenstates z/2, z/3 one has: 

z/^ = V2 cos 6*23 + z^3 siu 6*23 (6) 

= -1^2 sin ^23 + cos ^23 • 

If the masses m2 and ms are different, quantum mechanical time evolution of an initial 
state induces a non-zero transition probability to Vr- The survival probability for the muon 
neutrino is 

P{i^^. -> z/^) = 1 - sm\2d2^) ■ sm^ { j , (7) 

where L (in km) is the distance travelled by the neutrino, Ey (in GeV) its energy and 
Am^g = m| — m2 (in eV^). The three-flavour case is governed by two independent differences 
of mass squares and three mixing angles. The best-flt oscillation parameters derived from 
present data are [91j: 

I Amgi I = 2.40 X 10"^ eV^ ^ sin^2^23 ^ 1 

Am^i = 7.65 X 10~^ eV^ sin^ ^12 = 0.304 

sin^ < 0.056 (3(t) . (8) 
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If an atmospheric muon neutrino crosses the Earth along its full diameter, the first minimum 
of the survival probability occurs around 24GeV (see Fig. |6]), the second at 8 GeV. For 
shorter propagation paths through the Earth, the minima occur at correspondingly reduced 
energies and become unobservable when approaching the horizon. 



Muon neutrino survival probability 




Figure 6. The probability of muon neutrino 
disappearance as a function of neutrino energy 
El, and zenith angle for atmospheric neutrinos. 
A zenith of 180° refers to a vertically-upward 
going neutrino, i.e. propagation along the full 
Earth diameter. A zenith angle of 90° marks 
horizontal propagation with a baseline of a few 
100 km. The pattern at lower energies 
corresponds to "standard oscillations" (without 
taking into account matter effects), the 
patterns at higher energies would be expected 
for non-standard oscillations due to the 
violation of Lorentz invariance (see text). 



If their energy threshold is as low as about 10 GeV (like for IceCube's DeepCore, see 
Sect. [5]), underwater /ice detectors can approach the range of conventional oscillations at 
high energies [92j, complementary to underground detectors like Super-Kamiokande which 
are mostly sensitive to lower energies. 

For Galactic or extragalactic distance scales, kpc to Mpc and beyond, only average oscil- 
lation probabilities are observable since even tiny energy differences cause large oscillation 
phase shifts; in addition, the emission regions may well be larger than oscillation wave- 
lengths. It can be shown that the mixing parameters listed in eq. ^ transform a produced 
flavour ratio of z/g : t'/x : z^t = 1 : 2 : to an observed one of : : Ur = 1 : 1 : 1 |93t IM] . 

Matter effects and neutrino mass hierarchy: 

One of the remaining puzzles of neutrino physics is the hierarchy of the mass eigenstates 
1^1,1^2, 1^3- We know that m2 is larger than mi, but we do not know whether rris is the largest 
(normal hierarchy) or the smallest mass (inverted hierarchy). It is expected that eventually 
double-beta decay experiments may provide the answer (in case neutrinos are Majorana 
particles). Equally promising are long-baseline neutrino accelerator experiments, provided 
sin^ 2012, > 0.001. Also a 100 Megaton detector for neutrinos may give the answer if sin^ 2^13 
is not too much below 0.1 and the detection threshold is 5-10 GeV |92| l95| l96] . This is the 
energy range where matter- induced oscillations become important. 

For neutrinos below 15 GeV passing through the centre of the Earth, the MSW effect 
slightly enhances the z/^ oscillation probability and suppresses the oscillation probability. 
For the inverted hierarchy, in contrast, the oscillation probability for is enhanced and 
for z/^ suppressed. Water/ice neutrino telescopes cannot distinguish between neutrinos and 
antineutrinos (i.e. identify the charges of the secondary muons). However, in the relevant 
energy range the interaction cross sections of neutrinos and anti-neutrinos differ by about 
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a factor of two: cr(z/) ~ 2a{u). This translates into a difference in the number of observed 
muon events (about 5% for sin^ 2^13 = 0.1). Based on the statistical discrimination with 
respect to the number of events below 15 GeV, it may therefore be possible to distinguish 
normal from inverted hierarchy with DeepCore. Needless to say that firmly establishing 
such a small difference requires excellent knowledge of the detector systematics. Given the 
complicated optical properties of ice, it remains to be demonstrated whether systematic 
uncertainties for these low energies can be kept sufficiently small. In this respect, a densely 
instrumented water detector may be superior. 

Non-standard oscillations: 

Many quantum gravity (QG) models suggest that Lorentz invariance may be violated 
and that this effect may be seen in neutrino oscillations P7H^ . Conventional oscillations 
are due to different mass eigenstates. Violation of Lorentz invariance (VLI) may occur if 
in addition there are also different velocity eigenstates, each with its own limiting velocity, 
differing from the speed of light. The VLI magnitude depends on the velocity splitting 
Ac/c = (ci — C2)/c. Assuming maximum conventional and VLI mixing and the same phase, 
the muon survival probability can be written as 



where the first term comes from standard oscillations (see eq. ([7])) and the second from 
VLI mixing. Note that the two effects scale differently with energy, the one with 1/E, the 
other with E or, more generally, with i?" if Ac/c is replaced by a generalised VLI term 
A6 = Ac/c X E^-". 

Another possible consequence of QG is the evolution of pure neutrino states to mixed 
states via interaction with the foamy space-time itself (quantum decoherence). The resulting 
effect can be characterised by a set of parameters Di, . . . , which in the simplest case are 
equal and can vary with some integral power m of energy. 

Figure [6] shows the survival probability as a function of neutrino energy E^y and zenith 
angle for neutrinos created in the Earth atmosphere. The pattern at lower energies corre- 
sponds to "standard oscillations" (without taking into account matter effects), the pattern 
at higher energies would be expected for non-standard oscillations due to the violation of 
Lorentz invariance with Ac/c = 10'^''. 

2.3. Environmental and marine sciences 

The construction of neutrino telescopes in deep ice or deep water provides opportunities 
for long-term, real-time measurements in these environments that are unique and therefore 
of utmost interest to a variety of scientific disciplines beyond astroparticle physics, such as 
geology and geophysics, marine biology, oceanography or environmental sciences. Owing 
to the larger construction flexibility in water and to the richness of marine phenomena, 
corresponding activities are mostly concentrated on marine sites. 

Multidisciplinary observatories are associated to all current Mediterranean deep-sea neu- 
trino telescope projects (see |1UU[ IIUI) and |1U2] and references therein). For KM3NeT, a 




(9) 
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dedicated interface to such instrumentation is integral part of the planning |1U3| . and a close 
cooperation with the European Multidisciplinary Seafloor Observatory (EMSO) project |1U4] 
has been established. It is expected that deep-sea data will be taken by a series of instru- 
ments that are connected to the KM3NeT deep-sea cable network; in conjunction, also the 
neutrino telescope data themselves will be useful (e.g. light measurements for biolumines- 
cence studies, calibration data for acoustic and water current measurements). 
Issues to be addressed with these data are amongst others [103]: 

• Investigation of deep-sea phenomena such as internal waves and short time-scale oscil- 
lations in the water column, which are relevant for understanding the physical processes 
of the ocean and their effects on the distribution of suspended geological, chemical and 
biological materials; 

• Real-time tracking of bio-acoustic emissions or vertical migrations of organisms; 

• Long-term monitoring of the ocean margin environment around Europe; 

• Studies of oceanic processes and the land-ocean-atmospheric interactions. Geologi- 
cal records show that the ecosystem of the Mediterranean Sea amplifies climatologic 
variations, making it an ideal test bed for climate studies. 

In addition, contributions to hazard warning systems (in particular related to tsunamis) are 
conceivable, even though problems related to reliability standards, data access and usage 
rights need to be clarified. 

In the South Polar ice, the spectrum of multidisciplinary investigations is limited com- 
pared to Ocean detectors, but nevertheless manifold. Light propagation in deep ice is affected 
by remnant air bubbles at shallow depths |1U5] and by layered impurities from dust which 



are due to climatic effects or volcano eruptions |1U6] - see Fig. 19 in Sect. 4.3 Motivated 



by these findings, a project called Deeplce was proposed in 1999 to the US National Science 
Foundation. It was suggested to join efforts of astroparticle physicists, glaciologists, seis- 
mologists, geophysicists and biologists to investigate and exploit the features of glacial ice 
at or near the South Pole. The project was not funded, but many of its ideas have been 
pursued further and have led to relevant results, e.g. on the propagation of acoustic and 
radio waves in ice (see also Sect. [?]), on glaciology, climate research and volcanology. One 
legacy of Deeplce is an optical dust-logger, which reads out concentrations of dust particles 
and volcanic ash in the 2.5 km deep water-filled boreholes with a depth resolution of about 
1 mm [107] . These dust logs, together with similar results from Greenland ice, have revealed 
a relationship between abrupt climate changes and faint volcanic fallout layers and have 
permitted for reconstructions of dust and paleowind records of high quality. The detailed 
investigation of the IceCube dust logs have provided an unique characterisation of Polar 
deep ice and make the South Pole a leading candidate for a next US deep-ice core project. 

Moreover, the thermistors and microinclinometers installed in IceCube allow for studying 
the shear strain rate of a large volume of ice in three dimensions as a function of stress, 
impurity content and temperatures down to — 35°C. Last but not least, we mention an idea on 
studying biology under extreme conditions which also evolved from Deeplce: a quantitative 
model of how micron-sized microbes could live in liquid veins at triple junctions of ice grains 
and how they could be detected [108] . 
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3. Detection Principles 

The classical operation of neutrino telescopes underground, underwater and in deep ice 
is recording upward-going muons generated in charged current muon neutrino interactions. 
The upward signature guarantees the neutrino origin of the muon since no other known 
particle can traverse the Earth. Neutrino telescopes need to be situated at more than a kilo- 
metre depth in order to suppress downward-moving muons which may be misreconstructed 
as upward-moving ones (Fig. [?]). Apart from these, only one irreducible background to ex- 
traterrestric neutrinos remains: neutrinos generated by cosmic ray interactions in the Earth's 
atmosphere ("atmospheric neutrinos"). This background cannot be reduced by going deeper. 
On the other hand, it provides a standard calibration source and a reliable proof of principle. 




Figure 7. Sources of muons in deep underwater /ice 
detectors. Cosmic nuclei - protons (p), a particles (He), 
etc. - interact in the Earth atmosphere (light-coloured). 
Sufficiently energetic muons produced in these 
interactions ("atmospheric muons") can reach the 
detector (white box) from above. Upward-going muons 
must have been produced in neutrino interactions. 



Underwater /ice neutrino telescopes consist of an array of photomultiplier tubes (PMTs) 
housed in transparent pressure spheres which are spread over a large volume in oceans, 
lakes or glacial ice. The PMTs record arrival time and amplitude, sometimes even the full 
waveform, of Cherenkov light emitted by muons or particle cascades. In most designs the 
spheres are attached to strings which - in the case of water detectors - are moored at the 
ground and held vertically by buoys. The typical PMT spacing along a string is 10 -20 m, 
and the distance between adjacent strings 60 -150 m. The spacing is thus by orders of 
magnitude larger than in an underground detector like Super-Kamiokande. This allows for 
covering large volumes, but makes the detector practically blind with respect to events with 
energies below about 10 GeV. 

3.1. Neutrino interactions 

At neutrino energies E^, above some 10 GeV, as relevant for this review, charged-current 
(anti)neutrino-nucleon reactions are dominated by deep-inelastic scattering, ^i)\N — )■ i^X 
(charged current, CC) or VgN — )■ (neutral current, NC), with £ = e, fi, t. The leading- 

order differential cross section for the most important process, CC reactions, is 

favN^^^2Glm^ / y id{x,Q^) + il-y)M^,Q^) for z/ 

dxdy TT '[q^ + M^J '''\{l-yyu{x,Q') + d{x,Q') for z7 ^ ^ 
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where x and y are the Bjorken scahng variables (x being the fraction of the nucleon momen- 
tum carried by the struck quark, y the fraction of the incoming neutrino energy transferred to 
the hadronic system X), pa 2xyEyMi^ is the negative four-momentum transfer squared, 
Gi? is the Fermi constant and and My/ are the nucleon and the W mass, respectively (see 
|109] for more details). The parton distributions d{x,Q'^) and u{x,Q'^) represent the sums 
of all d-type and ti-type quark flavours, respectively; likewise d{x,Q'^) and u{x,Q^) are the 
corresponding antiquark distributions. These functions have been measured in fixed-target 
experiments and at HERA (see |110) and references therein). 

Integrating eq. (10) over x and y yields the total cross section (J^u'n^ij,^x relevant for 



neutrino telescope observations. For E,^ < lOTeV we have (Q^) ^ M^; in this regime the 
cross section is linear in E^, to a good approximation, with a value of a ~ 10~^^ cm^ at 
Ey = 1 TeV. For higher energies, the propagator term and the QCD-induced Q"^ dependence 
of the parton distributions result in a slower rise, roughly proportional to Ey''^. The cross 
sections resulting from a recent analysis [lllj for CC and NC reactions and their uncertainties 
are shown in Fig. [8j Note that the cross section is per nucleon and is different for neutrinos 
and antineutrinos; for isoscalar targets, the cross sections roughly have ratio v : V = 2 : 1, 
whereas they are approximately equal for proton targets. The uncertainties at large neutrino 
energies are mainly caused by the lack of experimental constraints on the parton distributions 
at X < 10"^ 
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Figure 8. Total cross sections for neutrino (left) and antineutrino (middle) scattering on isoscalar 
nucleons. Shown are results of two different analyses, both for charged current (CC) and neutral 
current (NC) scattering. In the right panel, the relative uncertainties in both analyses are indicated. 
Figures taken from |lTT] (CTW); shown are also results from [Il2j (GQRC) and [mj (CSS). 



The final state lepton follows the initial neutrino direction with a mean square root 
mismatch angle Q decreasing with the square root of the neutrino energy |21] : 



1.5° 



'Ey [TeV] 



In |114] . the parameterisation {^Q) ^ ^.1° j{Ey [TeV])'^'^ is given for the average angle. On the 
one hand, this principally allows for source tracing with CC muon neutrino reactions, but 
on the other hand represents an intrinsic kinematic limit to the ultimate angular resolution. 
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which is shghtly worse than for high-energy gamma-ray astronomy and orders of magnitude 
worse than for conventional astronomy. 

For a CC reaction of a muon neutrino with energy E,y, we define Py^^{Ey^E'^^^) to 
be the probabihty to produce a muon which reaches the detector with an energy exceed- 
ing the minimum detectable energy, i^^™™. This probability depends on the cross section 
dayN-^^x{.E^, E^)/dE^ and the effective muon range i?eff) which is defined as the range after 
which the muon energy has decreased to i?™'" 



P ( E R™'^- 



--IGeV) 



P 



N 



dE, 



da. 



{Ey, Efj 



dE„ 



Emm 



(12) 



with p being the density of the target material and M^- the nucleon mass. For water and 
1 GeV the following approximation holds : 



77' mm 



P, 



1.3 X 10-6 ■ (E^/lTeV)2-2 
1.3 X 10-6 ■ (E^/lTeV)°-s 



for Ey<l TeV 
for Ey > 1 TeV . 



(13) 



This implies that a neutrino telescope can detect a muon neutrino with 1 TeV energy with 
a probability of about 10^^ if the telescope is on the neutrino's path. 

For a neutrino flux ^y{Ey,i)) = [d'^N^/ldE^dtdAdQ)] ■ AQ arriving from a solid angle 
region AQ at zenith ■(9, the number of events recorded by a detector with area A within a 
time T is given by 



N,{Ef-,^) 
AT 



dE, ^,{E,, ^) ■ P.^,,{E,, Ef^) ■ exp 



(^'-v-N~>t,^x{E,) ■ Z{d) 



M, 



N 



■ (14) 



77min 



Here Z{d) is the matter column density in the Earth crossed by the neutrino, in g/cm^. For 
sub- TeV energies, absorption in the Earth is negligible and the exponential term is close to 
unity; for larger neutrino energies, absorption becomes relevant (see Fig. [9]). 
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Figure 9. Transmission probabihty 
through the Earth for neutrinos of 
different energies, as a function of the 
zenith angle. For each energy, the upper 
line is for CC interactions and the lower 
one for CC and NC interactions. Figure 
taken from )115| . 



The fraction of the neutrino energy carried by the final-state lepton is given by 1 — y. For 



lower energies, where {Q ) <^ M^, integration of eq. (10) yields 1/4 < {y) < 1/2, depending 
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on the relative contribution of quarks and antiquarks. For larger Ei^, the dependence of the 
propagator term on ^ 2xyM]<^Ei, reduces the mean y (see Fig. 10). The muon in z/^ 
CC reactions on average thus carries the major fraction of the neutrino energy, whereas the 
hadronic cascade has lower energy. 
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Figure 10. Average y as a function of 
neutrino energy, for CC (solid lines) and 
NC (dashed) reactions. Figure taken 
from |116) . 



3.2. Cherenkov light 

Charged particles moving faster than the speed of light in a medium with index of 
refraction n, i.e. v > c/n, emit Cherenkov light under the characteristic angle 

cos 6'c = "5^, (15) 
p ■ n 

where (i = v/c and n depends on the frequency v of the emitted photons, n = n{v). The 
spectral distribution of Cherenkov photons per path length of an emitting particle with 
charge ±2;e is given by 

dA^ ( 1 



dxdA A2 /32.n2j ^^^^ 

with a being the fine structure constant. The total amount of released energy per particle 
path length is obtained by multiplying eq. (16) with = hu and integrating over u: 



l3-n(u)>l 

For A = 550 nm the index of refraction for water is n ~ 1.33, yielding about 400 eV/ cm, 
or approximately 200 Cherenkov photons per cm in the transparency window of water, i.e. 
for wavelengths 400 nm < A < 700 nm. The average Cherenkov angle of these photons is 
Be ^ 43°. 

3. 3. Light propagation and detection 

The propagation of light in water is governed by absorption and scattering, which both 
depend on the wavelength A. In the first case the photon is lost, in the second case it changes 
its direction. Scattering effectively delays the propagation of photons between the points of 
emission and detection. The parameters generally chosen to describe these phenomena are: 
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1. The absorption length La{X) - or the absorption coefficient a(A) = l/L^ - describes 
the exponential decrease of the number of non-absorbed photons as a function of 
photon path length r, N = Nq ■ exp(— r/La). 

2. The scattering length Lfc(A) and scattering coefficient 6(A) = 1/Lfy, defined in analogy 
to La{X) and a(A). 

3. The scattering function x(6', A), i.e. the distribution in scattering angle 6. 

4. Instead of the "geometrical" scattering length Lh{X) one frequently uses the effective 
scattering length Lgg = Liy/{1 — (cos 6')), where (cos 6') is the mean cosine of the 
scattering angle. Lgg "normalises" scattering lengths for different distributions A) 
of the scattering angle to the extreme case {cos 6) = 0, i.e. Les can be interpreted as 
isotropisation length. For {cos 6) ~ 0.8-0.95, as for all media considered here, photon 
delay effects in media with the same Les are approximately the same. 

In Table [T| typical parameter values for Lake Baikal, oceans and the Antarctic ice are 
summarised (see |117| for references). All values are given for the wavelength where they 
are maximal. 



Site 


La [m] 


^eff [ni] 


Lake Baikal, 1 km depth 


22 


150-400 
(seasonal variations) 


Ocean, > L5km depth 


40-70 (depends 
on site and season) 


200-400 (depends 
on site and season) 


Polar ice, 1.5 -2.0 km depth 


^ 95 (average) 


^ 20 (average) 


Polar ice, 2.2 -2.5 km depth 


~ 150 (average) 


~ 40 (average) 



Table 1. Absorption length and effective scattering length for deep lake and ocean water and for 
deep ice. 



Scattering and absorption in water and ice are measured with artificial light sources. The 
scattering coefficient, to a large part due to particulate matter, changes with wavelength 
less pronouncedly than the absorption coefficient (see Fig. 20 for ice). In water, the depth 
dependence over the vertical dimensions of a neutrino telescope is small, but parameters 
may vary with time, due to transient water inflows loaded with bio-matter or dust, or due 
to seasonal changes in water parameters. They must therefore be permanently monitored. 
In glacial ice at the South Pole, the situation is different. The parameters are constant in 
time but strongly depend on depth (see below). 

Strong absorption leads to reduced photon collection, strong scattering deteriorates the 
time information which is essential for the precise reconstruction of tracks and showers (see 
below) . 

Efficient and accurate detection and reconstruction of neutrino events requires light de- 
tection with a sensitivity at the single-photon level and a measurement of the arrival time 
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with nanosecond precision. Furthermore, the integrated photo-sensitive area, multiphed 
with the single-photon detection efficiency, is the prime parameter determining the sensitiv- 
ity of the neutrino telescope and is to be maximised. Currently, PMTs |118l 1119) are the 
only devices matching these constraints at an affordable price. Relevant parameters are the 
photocathode size (i.e. the diameter of the tube); the quantum efficiency (probability that 
an incoming photon causes emission of a photo-electron) and its spectral dependence; the 
collection efficiency (probability that a photo-electron causes a measurable signal); the tran- 
sit time spread (jitter of the delay between photon interaction and output signal); the gain 
(electric amplification, i.e. average number of electrons in the signal of one photo-electron). 

• ooooooooo^oooo 

PMTs 




Cherenkov cone spherical Cherenkov front 

• ooooooooooooo 



Figure 11. Detection principles for muon tracks (left) and cascades (right) in underwater /ice 
detectors. Note that the Cherenkov light emission by cascades is peaked at the Cherenkov angle 6c 
with respect to the cascade axis but has a wide distribution covering the full solid angle. 



3.4- Detection of muon tracks and cascades 



Figure 11 sketches the two basic detection modes of underwater /ice neutrino telescopes. 
In the following, denotes the target nucleon and X the hadronic final state. CC muon 
neutrino interactions, "u'^N — )■ fi^X produce a muon track (left; an IceCube event display 



of a 10 TeV muon is shown in Fig. 12), whereas other neutrino reaction types cause hadronic 
and/or electromagnetic cascades (right). This is in particular true for NC reactions (VW — )■ 
hadronic cascade) or CC reactions of electron neutrinos i^y' — )■ e^X, overlapping 
hadronic and electromagnetic cascades). CC tau neutrino interactions, V^-N — )■ r^X, can 
have either signature, depending on the r decay mode. In most astrophysical models, neu- 
trinos are produced with a flavour ratio i/g • ^'u ^ '^r ~ 1 • 2 : 0. Over cosmic distances. 



oscillations turn this ratio to i/g : ^^/^ : J^r ~ 1 : 1 : 1 (cf. Sect. 2.2.3), which means that almost 
2/3 of the charged current interactions initiate cascades only. 

3.4- i- Muon tracks 

The most precise measurement of the neutrino direction is achieved for charged current 
muon (anti)neutrino interactions; this channel is therefore central to all investigations of 
astrophysical neutrino sources. The muon direction, and also its energy, are reconstructed 
from the Cherenkov cone accompanying upward-going muons. The upward signature guar- 
antees the neutrino origin of the muon since no other particle can cross the Earth. The 
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Figure 12. Event displays of a muon (left) and a cascade (right) event recorded in IceCube. The 
coloured circles indicate the photomultipliers hit; their sizes correspond to the numbers of photons 
detected. The colour code visualises the hit times, with red hits earliest and violet ones latest. 
The reconstructed zenith angle and energy of the muon are 117° and 20TeV, respectively; the 
reconstructed cascade energy is 175 TeV. 



effective target volume considerably exceeds the actual detector volume due to the large 
range of muons (about 0.96 x lO^g/cm^ (1.04 x lO^g/cm^) at 300 GeV and 14 x lO^g/cm^ 
(12 X lO^g/cm^) at 100 TeV in water/ice (standard rock), logarithmically rising towards 
larger energies |120) : lO^g/cm^ corresponds to 1km in pure water). 

The muon loses energy mainly via ionisation, pair production, bremsstrahlung and photo- 
nuclear reactions. The energy loss can be parameterised by [21^ il20j 

-'^ = a + b-E,, (18) 

where a and b are weakly energy dependent. For water, the ionisation loss is roughly given 
by a = 2MeV/cm. The energy loss from pair production, bremsstrahlung and photo- 
nuclear reactions becomes dominant at energies beyond 1 TeV, where it is approximately 
linear in E^; for water and = 100 TeV it is given by 6tot = ^pair + ^brems + ^phonuc = 
(1.7 + 1.2 + 0.6) X 10"^gcm"2 = 3 5 x lO^^gcm'^ jT^O]- Note that all energy loss mech- 
anisms producing relativistic charged particles - in particular from pair production and 
conversion of bremsstrahlung photons - also lead to additional Cherenkov light. 

For the reconstruction of a muon track, its Cherenkov photons need to be measured at 
least at five positions. In reality, because of timing inaccuracies due to scattering (mostly in 
ice) and to light background from radioactive decays (deep sea) and bioluminescence (lakes 
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and deep reliable reconstruction requires typically 10 or more PMT signals. 

Underwater/ice telescopes are optimised for the detection of muon tracks at energies of 
the order TeV or above, for the following reasons: 

1. The flux of neutrinos from cosmic accelerators is expected to be harder (generic 
shape E~'^) than that of atmospheric neutrinos {E~^''^), yielding a better signal-to- 
background ratio at higher energies. 

2. Neutrino cross section and muon range increase with energy. The larger the muon 
range, the larger the effective detector volume (even though this argument must be 
used with care since a muon that has lost most of its energy when reaching the detector 
may be discarded due to energy-dependent selection cuts). 



3. The mean angle between muon and neutrino decreases with energy (see eq. (11)), 
resulting in better source tracing and signal-to-background ratio at high energy. 

4. For energies above a TeV, the light emission increases linearly with Efj_, providing an 
estimate for the muon energy with an accuracy of cr(logi?^) ~ 0.3. By unfolding pro- 
cedures, a muon energy spectrum can be translated into a neutrino energy spectrum. 

Muons which have been generated by cosmic ray interactions in the atmosphere above 
the detector and punch through the water or ice down to the detector outnumber neutrino- 
induced upward-moving muons by several orders of magnitude (about 10^ at 1 km depth 
and 10"^ at 4 km depth) and have to be removed by careful up/down discrimination. 

At energies above a few hundred TeV, where the Earth becomes opaque even to neutrinos, 
neutrino-generated muons arrive preferentially from directions close to the horizon, at EeV 
energies essentially only from the upper hemisphere (Fig.|9]). The high energy deposition of 
muons from PeV-EeV extraterrestrial neutrinos provides a handle to distinguish them from 
downward-going atmospheric muons (that have a spectrum decreasing much more steeply 
with energy). 

3.4-2. Cascades 

Charged current interactions of electron and, in most cases, tau neutrinos and all neutral 
current interactions do not lead to high energy muons but to electromagnetic or hadronic 
cascades. Their length increases only with the logarithm of the cascade energy and is of the 
order 5 - 20 m. Cascade events therefore typically need to be contained in the instrumented 
volume of the detector to be reconstructible ("contained events"). 

With their small length and a diameter of the order of 10 - 20 cm, cascades may be 
considered as quasi point-like compared to the spacing of the PMTs. The effective volume 
for cascade detection is close to the instrumented volume of the detector. For first-generation 
neutrino telescopes it is therefore much smaller than that for muon detection. However, for 
kilometre-scale detectors and moderate energies both can be of the same order of magnitude. 
The total amount of light is proportional to the cascade energy and thus provides a direct 
measurement of E for contained cascades, in contrast to muons, for which only dE/dx 
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can be measured. Therefore, in charged current z/g interactions, the neutrino energy can 
be determined with an accuracy of about 30% (note, however, that this is not true for NC 
events; also, in CC events the Ur from the decay of the final-state r carries away a sizeable 
fraction of "invisible" energy). Whereas the energy resolution is much better than for muons, 
the directional accuracy is worse since the lever arm for fitting the direction is small (for 
actual values, see below). The irreducible background from atmospheric electron neutrinos 
is significantly smaller than that of atmospheric muon neutrinos for the muon channel. All 
this makes the cascade channel particularly interesting for searches for a diffuse, high-energy 
excess of extraterrestrial over atmospheric neutrinos. 

3.5. Effective area and sensitivity 

The minimum energy of a muon to be recorded depends on the impact point and the 
impact angle, resulting in a smearing for -E™™. Even muons outside the geometrical area 
may fulfil the trigger condition. On the other hand, background suppression usually requires 
tight cuts during the analysis, so that muons having triggered the detector might be rejected 
at a later stage. Therefore, the geometrical area A is replaced by the effective muon area 
Agjj which is calculated from Monte Carlo simulations as follows: N^^n muons are generated 
perpendicular to an area A^en significantly larger than the detector cross section. With 
A^passed muous acccptcd after all cuts, the effective area is then given by 



= ^gen ■ . (19) 

'generated 

The effective area depends on 6 and 0, i.e. the muon direction with respect to the detector. 
Due to the increasing range and light emission, increases with muon energy. Note that 
A'^Q also depends on the distance between the point where the muons are generated and the 
detector. 

The concept of an effective area can also be applied to neutrinos, taking into account the 
neutrino nucleon cross section, the muon range and the absorption in the Earth according 



to Sect. |3.1[ Figure [13] shows the effective neutrino area of IceCube as a function of energy 
and zenith angle. The strong increase with energy is due to the increasing neutrino cross 
section and to the increasing muon range. The decrease at high energies for large zenith 
angles is due to neutrino absorption in the Earth. 

Due to the small neutrino interaction cross section, the effective neutrino area is many 
orders of magnitude smaller than the muon effective area. Even cubic-kilometre neutrino 
telescopes reach only neutrino areas between a few square metres and a few hundred square 
metres, depending on energy. This has to be compared to several ten thousand square metres 
typical for air Cherenkov telescopes which detect gamma-ray initiated air showers. A ratio 
1:1000 (10 m^ : 10000 m^) may appear desperately small. However, one has to take into 
account that Cherenkov gamma telescopes can only observe one source at a time, and that 
their observations are restricted to Moon-less, clear nights. Neutrino telescopes observe a full 
hemisphere, 24 hours per day. Therefore, cubic-kilometre detectors reach a flux sensitivity 
similar to that reached by first-generation Cherenkov gamma telescopes like Whipple and 
HEGRA [1221 [123! for TeV gamma rays, i.e. $(>1 TeV) ^ lO^^^ ^,^^2 g-i 
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zenith range (90°, 120°) 
zenith range (120°, 150°) 
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log^^ E„ (GeV) 



Figure 13. Effective neutrino area 
of the 40-string configuration of 
IceCube as a function of energy, 
averaged over different zenitli angle 
ranges. At energies beyond 
10^ GeV, neutrino absorption n tlie 
Earth becomes important. Figure 
taken from |121j . 



3. 6. Reconstruction 

In this section, some relevant aspects of event reconstruction are discussed for the case 
of muon tracks |124| - [126] . For cascades, see |127H129] . The reconstruction procedure for a 
muon track consists of several consecutive steps which are typically: 

1. Rejection of noise hits; 

2. Simple pre-fit procedures providing a first-guess estimate for the following iterative 
maximum-likelihood reconstruction; 

3. Maximum-likelihood reconstruction; 

4. Quality cuts in order to reduce background contaminations and to enrich the sample 
with signal events. This step is strongly dependent on details of the actual analysis - 
diffuse fluxes at high energies, searches for steady point sources, searches for transient 
sources etc. 



An infinitely long muon track can be described by an arbitrary point tq on the track 



which is passed by the muon at time to? with a direction p and energy Eq (see Fig. 14) 



Since the muon in very good approximation propagates with the vacuum speed of light, c. 



photons emitted under the Cherenkov angle 6^ (see eq. (15)) and propagating on a straight 



path ("direct photons") are expected to arrive at PMT i located at r*j at a time 

p-ifi — roi + d- tan 6c 
tgeo = to -\ , (20) 

where d is the closest distance between PMT i and the track. The time residual tres is given 
by the difference between the measured hit time thit and the hit time expected for a direct 
photon, tgeo- 

^res ^hit ^geo • (^-^) 
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Y, p=1/n 



A 
P 



(ro,to) 



V i ' i / Figure 14. Geometry of muon reconstruction: The 
muon, propagating with /3 = 1 in very good 
approximation, radiates a Cherenkov photon under 
the angle 6c (with cos Be = ^/n) that propagates with 
/3 = 1/n on a straight hne and hits the photon sensor 
at position rj. Under these assumptions the hit time 
can be calculated to be ti = tgeo (see eq. (20)). The 
time difference ti — to thus relates the muon track 
parameters (position vq, direction p) to the position 
of the photon hit. 



Schematic distributions for time residuals are given in Fig. 15 An unavoidable symmetric 
contribution around = in the range of a nanosecond comes from the PMT/electronics 
time jitter, ctj. An admixture of noise hits to the true hits from a muon track adds a flat 



pedestal contribution like shown in top right of Fig. 15 Electromagnetic cascades along the 
track, initiated by bremsstrahlung and pair production, lead to a tail towards larger (and 
only larger) time residuals (bottom left). Scattering of photons can lead to an even stronger 
delay of the arrival time (bottom right). These residuals must be properly implemented 
in the probability density function for the arrival times used in the maximum-likelihood 
procedure. 
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Figure 15. Schematic distributions of arrival times 
for different cases (see text). 
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The simplest likelihood function is based exclusively on the measured arrival times. It is 
the product of all A^hit probability density functions Pi to observe, for a given value of track 
parameters {a}, photons at times U at the location of the PMTs hit: 

A^hit 

Ltua.^\[p{tr..,i\{(^}) ■ (22) 
1=1 

More complicated likelihood functions include the probability of hit PMTs to be hit and 
of non-hit PMTs not to be hit, or of the respective amplitudes. Instead of referring only to 
the arrival time of the first photon for a given track hypothesis and the amplitude for a given 
energy hypothesis, one may also refer to the full waveform from multiple photons hitting the 
PMT. For efficient background suppression, the likehhood may also incorporate information 
about the zenith angular dependence of background and signal (Bayesian probability). The 
reconstruction procedure finds the best track hypothesis by maximising the likelihood. 
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4. First-Generation Neutrino Telescopes 



An explorative phase of more than two decades was required to solve the technical 
problems of instrumenting target masses of the order of 10 Megatons in deep ice, fresh or sea 
water and to overcome the difficulties which these hostile environments pose. Nevertheless, 
neutrino telescopes of this size have been constructed and successfully operated in all three 
media. In this Section, the technical setup of these installations is discussed, also with 
reference to the history of major steps leading to their construction. 

4.1. DUMAND 

The history of underwater neutrino telescopes starts with a project which eventually 
failed but left an incredibly rich legacy of ideas and technical principles: The DUMAND 
project. DUMAND stands for Deep Underwater Muon and Neutrino Detector. Its early 
history is excellently covered in a "Personal history of the DUMAND project" by A. Roberts 
|130] . At the 1973 International Cosmic Ray Conference (ICRC), a small group of physicists 
conceived a deep-water detector to clarify puzzles in muon depth-intensity curves. The 
puzzles faded away, but it was obvious that such a detector could also work for neutrinos. 
The year 1975 saw the first of a - meanwhile legendary - series of DUMAND Workshops. 
Soon, the decision was taken to deploy the detector 30 km off the coast of Big Island, Hawaii, 
at a depth of 4.8 km. 



The original idea to construct a cubic-kilometre detector (1978) with more than 20 000 
large-size photomultipliers (see Fig. [I6| was quickly abandoned due to technical and financial 
reasons. A half-sized configuration (1980) met the same fate, as did a much smaller array 
with 756 phototubes (1982). The latter design was comparable in size to the meanwhile 



1978: 1.26 

22,698 OMs 



km3 1988: 0.002 km3 
216 OMs 




Figure 16. The 

originally conceived 
DUMAND 
cubic-kilometre 
detector and the 1988 
plan for the 
first-generation 
underwater neutrino 
telescope DUMAND-II. 
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decommissioned AMANDA detector at the South Pole (see Sect. 4.3) and the operating 
ANTARES telescope in the Mediterranean Sea, close to Toulon (see Sect. 4.4). What finally 
emerged as a technical project was a 216-phototube version, dubbed DUMAND-II or "The 
Octagon" (eight strings at the corners of an octagon and one in the centre), 100 m in diameter 
and 230 m in height [131j. 

Signals on at least three strings are required for full spatial reconstruction of a muon 
trajectory. Consequently, the sequential deployment was planned to start with TRIAD, a 
sub-array of three full strings. Unfortunately, pressure housings of the first string developed 
leaks during deployment in 1993 and soon the communication to shore failed. The "3-string 
race" between DUMAND, the Baikal project and AMANDA was eventually won by the 
Baikal collaboration in 1993, whereas the DUMAND project was terminated in 1995. 

Ironically, the most sustaining physics result from the DUMAND project has been ob- 
tained with a 7-phototube test string deployed for only some hours from a ship |132) . It 
gives the measured muon intensity as a function of depth and thus returned to the initial 
idea of the 1973 ICRC. 

4-2. The Baikal neutrino telescope 

The Baikal Neutrino Telescope is installed in the Southern part of Lake Baikal |1331I134] . 
The distance to shore is 3.6 km, the depth of the lake at this location is 1366 m, the active 
part of the detector is located at a depth of about 1.1 km. 

First site surveys started in 1980. In 1984, a first stationary string was deployed and 
muons recorded |135| . This was followed by another stationary string in 1986 which was 
optimised for the detection of magnetic monopoles catalysing proton decay |136[ 1137] . 

A new period began with the development of the QUASAR photodetector (see below), 
which replaced the former 15 cm fiat photomultipliers, and the design of the array NT200 



see Fig. 17, left). The BAIKAL collaboration was not only the first, in 1993, to deploy three 



strings (as necessary for full spatial reconstruction of muon trajectories), but also reported 



the first atmospheric neutrino detected underwater (Fig. 17, right). 

The central part of the Baikal configuration is NT200, an array of 192 optical modules 
which was completed in April 1998 and has been taking data since then. The optical modules 
are carried by eight strings which are attached to an umbrella-like frame. The strings are 
anchored by weights at the lake fioor and held in a vertical position by buoys at various 
depths. The configuration spans 72 m in height and 43 m in diameter. The detector is 
deployed (or hauled up for maintenance) within a period of about 6 weeks in February to 
April, when the lake is covered with a thick ice layer providing a stable working platform. 
It is connected to shore by several copper cables on the lake fioor which allow for operation 
over the full year. 

The optical modules are glass spheres equipped with QUASAR-370 phototubes; they 
are grouped pair-wise along a string. In order to suppress accidental hits from dark noise 
(about 30 kHz) and bioluminescence (typically 50 kHz but seasonally raising up to hundreds 
of kHz), the two photomultipliers of a pair are switched in coincidence, defining a channel, 
with typically only 0.1 kHz noise rate. The basic cell of NT200 consists of a svjaska (Russian 
for "bundle"), comprising two optical module pairs and an electronics module for time and 
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Figure 17. Left: The 
Baikal Neutrino Telescope 
NT200. Right: One of the 
first upward moving muons 
from a neutrino interaction 
recorded with the 4-string 
stage of the detector in 
1996 [T38J. The Cherenkov 
light from the muon is 
recorded by 19 channels. 



amplitude conversion and slow control functions (Fig.[T7| left). A majority trigger is formed 
if > 4 channels are fired within a time window of 500 ns (this is about twice the time a 
relativistic particle needs to cross the NT200 array). Trigger and inter-string synchronisation 
electronics are housed in an array electronics module at the top of the umbrella frame. 
This module is less than 100 m away from the optical modules, allowing for nanosecond 
synchronisation over copper cable. 
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Figure 18. Left: The QUASAR-370 phototube. Right: A fuh Baikal optical module. 



Figure 18 shows the phototube and the full optical module [139J. The QUASAR-370 is a 
hybrid device. Photoelectrons from a large hemispherical cathode (K2CsSb) with > 27r view- 
ing angle are accelerated by 25 kV to a fast, high-gain scintillator which is placed near the 
centre of the glass bulb. The light from the scintillator is read out by a small conventional 
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photomultiplier (type UGON). One photo-electron from the hemispherical photocathode 
yields typically 20 photoelectrons in the conventional photomultiplier. This high multipli- 
cation factor results in an excellent single electron resolution of 70%. Furthermore, the 
QUASAR-370 is characterised by a small time jitter (2 ns) and a small sensitivity to the 
Earth's magnetic field. 

The small spacing of modules leads to a comparably low energy threshold of about 
15 GeV for muon detection. About 400 upward muon events were collected over 5 years. 
Still, NT200 could compete with the much larger AMANDA for a while by searching for 
high energy cascades below NT200, surveying a volume about ten times as large as NT200 
itself |140| . In order to improve pattern recognition for these studies, NT200 was fenced 
in 2005-2007 with three sparsely instrumented outer strings (6 optical module pairs per 
string). This configuration is named NT200+ |141] . 

4.3. AMANDA 

AMANDA (Antarctic Muon And Neutrino Detection Array) used the 3 km thick ice layer 
at the South Pole as target and detection medium |142l I143j . AMANDA was located some 
hundred metres away from the Amundsen- Scott station. Holes of 60 cm diameter are drilled 
with pressurised hot water; strings with optical modules are deployed in the molten water 
which subsequently refreezes. Installation operations at the South Pole are performed in 
the Antarctic summer, November to February. For the rest of the time, two operators (of a 
winter-over crew of 25-40 persons in total) maintain the detector, connected to the outside 
world via satellite communication. 
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Figure 19. The AMANDA configuration. Three of 
the 19 strings have been sparsely equipped towards 
larger and smaller depth in order to explore ice 
properties, one string got stuck during deployment at 
too shallow depth and was not used in analyses. The 
Eifi'el tower is shown to scale for size comparison. 
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Figure 19 shows the AMANDA configuration. A first test array with 80 optical modules 
at four strings (not shown) was deployed in 1993/94, at depths between 800 and 1000 m 
[105] . The effective scattering length Les was found to be extremely small, between 40 cm 
at 830 m depth and 80 cm at 970 m. The scattering is due to air bubbles trapped in the 
ice and makes track reconstruction impossible. The tendency of the scattering to decrease 
with depth, as well as results from ice core analyses at other places in Antarctica, suggested 
that the bubbles should disappear below 1300 m. This expectation was confirmed with a 
second 4-string array which was deployed in 1995/96. The remaining scattering, averaged 
over 1500-2000 m depth, corresponds to Les ~ 20 m and is assumed to be due to dust. This 
is still considerably worse than for water but sufficient for track reconstruction |106l I124| . 
The array was upgraded stepwise until January 2000 and eventually comprised 19 strings 
with a total of 677 optical modules, most of them at depths between 1500 and 2000 m. 
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Figure 20. Absorption coefficient (left) and scattering coefficient (right) in the South Polar ice as 
functions of depth and wavelength. 



In Fig. 20, absorption and scattering coefficients are shown as functions of depth and 



wavelength |106] . The variations with depth are due to bubbles at shallow depth leading 
to very strong scattering and, at larger depths, to dust and other material transported to 
Antarctica during varying climate epochs. The quality of the ice improves substantially 
below a major dust layer at a depth of about 2000 -2100 m, with a scattering length about 
twice as large as for the region above 2000 m. The depth dependence of the optical properties 
complicates the analysis of the experimental data. Furthermore, the large delays in photon 
propagation due to the strong scattering cause a worse angular resolution of deep-ice detec- 
tors compared to water (see below). On the other hand, the large absorption length, with a 
cut-off below 300 nm instead of 350-400 nm in water, results in better photon collection. 

The relatively short distance between optical modules and surface electronics allowed 
for transporting the analogue signals of the photomultipliers to surface over 2 km of cable 
instead of digitising them in situ. This requires a large output signal of the photomultiplier. 
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a specification met by the 8-inch R5912-2 from Hamamatsu with 14 dynodes and a gain 
of 10^. The first ten strings used copper cables for both high-voltage supply and signal 
transmission; for the last 9 strings the anode signal was used to drive a LED and the light 
signal was transmitted to surface via optical fibre. The time resolution achieved was 5-7 ns; 
given the strong smearing of photon arrival times due to light scattering in ice, this jitter 
appeared to be acceptable. An event was defined by a majority trigger formed in the surface 
counting house, requesting > 8 hits within a sliding window of 2 \xs. 

Time calibration of the AMANDA array was performed with a YAG laser at surface, 
sending short pulses via optical fibres to each optical module. This laser system was also used 
to measure the delay of optical pulses propagating between strings and so to determine the 
optical ice properties as well as the inter-string distances. A nitrogen laser (337 nm), halogen 
lamps (350 and 380 nm) and LED beacons (450 nm) located in deep ice yielded further 
information about the ice properties across a large range of wavelengths. The measured 
time delays were fitted and the resulting parameterisations implemented in the probability 



density functions for the reconstruction procedure (see Sect. 3.6). 

A big advantage compared to underwater detectors is the small photomultiplier noise 
rate, about 1 kHz in an 8-inch tube, compared to 20-40 kHz due to K'^° decays and biolumi- 
nescence in lakes and oceans. The contamination of hit patterns from particle interactions 
with noise hits is thus small and makes hit selection much easier than in water. 

The angular resolution of AMANDA for muon tracks is 2° -2.5°, with a lower energy 
threshold around 50GeV. Although better than for Lake Baikal (3° -4°), it is much worse 
than for ANTARES (< 0.5°, see below). This is the result of the strong light scattering 
which deteriorates the original information contained in the Cherenkov cone. The effect is 
even worse for cascades, where the angular resolution achieved with present algorithms is 
only 25° (compared to 5° -8° in water [128j ). 

4.4. ANTARES 



Starting with the decline of the DUMAND project (see Sect. 4.1), the participating Eu 



ropean groups started to explore options for a deep-see neutrino telescope in the Mediter- 



ranean Sea. The NESTOR (see Sect. |4^, ANTARES and NEMO (see Sect. |43]) projects 
were initiated, in the temporal order indicated. The NESTOR and NEMO Collaborations 
have performed technical R&D work and deployed prototype installations; they are now, to- 
gether with the ANTARES collaboration, pursuing further work in the KM3NeT framework 



(see Sect. 5.2). ANTARES has succeeded in installing and is operating the first working 
deep-sea neutrino telescope, thus providing the proof of feasibility for such devices. 

The ANTARES (Astronomy with a Neutrino Telescope and Abyss environmental RE- 
Search) |144) proposal |114] was presented in 1999. It was based on the operation of a 
demonstrator string |145l 1146) as well as on the results of extensive site exploration cam- 
paigns in the region off Toulon at the French Mediterranean coast, indicating that the optical 
background [147J as well as sedimentation and biofouling |148| are acceptable at that site. 
The ANTARES design encompasses 12 strings, each carrying 25 "storeys" equipped with 
three optical modules, an electronics container and calibration devices where necessary. A 
further string, the "instrumentation line", carries devices for environmental monitoring. The 
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strings are anchored to the sea floor with a dead weight and kept upright by a submerged 
buoy at their top. The cable connecting the storeys encompasses copper leads and optical 
fibres and at the same time has to sustain the mechanical tension; this combination was 
found to be conceptually challenging. At the sea floor, electro-optical cables connect the 
strings to a junction box, from which the main cable goes to shore. The strings are deployed 
from the surface and subsequently connected to the junction box employing a submersible. 
The inter-string distances are about 70 m, the vertical distance between adjacent storeys 
is 14.5 m. The depth at the ANTARES site is 2475 m. The schematic setup is shown in 



Fig. 21 , a detailed technical description can be found in |149 
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Figure 21. Schematic 
of the ANTARES 
detector. Indicated are 
the 12 strings and the 
instrumentation hne in 
its 2007 configuration 
(IL07). Shown as an 
inset is the photograph 
of a storey carrying 3 
photomultipHers . 



The ANTARES construction started in 2002 with the deployment of the main cable and 
the junction box. In 2002/2003, a preproduction string was deployed and operated for some 
months. Several technical problems were identifled that required further studies, design 
modiflcations and the operation of a mechanical test string |150] . The detector in its flnal 
conflguration was eventually installed in 2006-2008 and has been operational since then, 
with a break of a few months in 2009 due to a failure of the main cable that required repair. 

The ANTARES optical module consists of a 17-inch glass sphere housing a hemi- 
spherical 10-inch photomultiplier (Hamamatsu R7081-20), connected to the glass surface by 
optical gel, and its high-voltage base. A mu-metal cage shields the photomultiplier against 



the Earth magnetic fleld (see Fig. 22). A single cable with copper leads is used for voltage 



supply and for transporting the analogue photomultiplier signals to the digitisation elec- 
tronics in the electronics container of the storey; dispersion effects are negligible due to 
the short cable length of about a metre. The optical modules are oriented such that the 
photomultipliers look downward at an angle of 45°. 



36 



Figure 22. Photograph of an ANTARES optical 
module. Its upper hemisphere is painted black to 
reduce background light. 



The photomultipliers provide an intrinsic timing resolution (transit time spread) of about 
c"tts = 1-3 ns [152J, thus allowing for precise measurements of the arrival times of Cherenkov 
photons. Owing also to the minute light scattering in deep-sea water |153| and accurate 
timing and optical module position calibration (see below), muon trajectories can be recon- 
structed with a precision of about 0.2° for muon energies exceeding 1 TeV |154) . Degradations 
are expected for down-going muons (where the direct Cherenkov light is partly outside the 
photomultiplier acceptance and the fraction of scattered light is therefore increased) and 
due to background light from decays and bioluminescence. The single photon rate from 
decays is about 30 kHz per photomultiplier, the bioluminescence rate exhibits a slowly 
varying baseline of typically a few 10 kHz to a few 100 kHz, and second-scale "bursts" in the 
MHz region. Efficient data taking is possible up to baseline rates of 200 kHz. 

The ANTARES data acquisition |155j follows the "all-data-to-shore" concept, i.e. all 
photomultiplier signals exceeding an adjustable threshold (default 0.3 photo-electrons) are 
read out, thus avoiding any inter-storey or inter-string trigger processing off-shore. The 
signals are digitised and time-stamped by the storey electronics and then sent to shore via 
optical fibres. The resulting data rate is dominated by background hits and exceeds by far 
the capacity of data mass storage. Data filter software, running on an on-shore farm of PCs, 
selects event candidates based on multi-hit coincidences and/or multi-photo-electron hits in 
single photomultipliers. The selected data are stored for offline analysis. 

The two main calibration tasks are the synchronisation of the signal time measurements 
at the individual photomultipliers and the position and orientation monitoring of the optical 
modules that move with their strings in the sea current. For the timing calibration [156J, 
light signal running times through the optical fibres from shore to each storey and back 
are used, as well as pulsed light emissions by laser and LED beacons |157) . An accuracy 
of about 0.5 ns is achieved, complying with the requirements. The position and orientation 
calibration |158t 1159] uses data from compasses and tiltmeters on each storey and from an 
acoustic system measuring running times of acoustic pulses between transmitters on the sea 
fioor and receivers (hydrophones) on 5 storeys per string. Also here, the precision of roughly 
10 cm, translating into a timing uncertainty of 0.5 ns, is within specifications. 

ANTARES is equipped with an acoustic detection system, AMADEUS |160] . for feasibil- 



ity studies towards acoustic neutrino detection. For a more detailed discussion see Sect. 7.3 
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4.5. NEMO 

The NEMO (NEutrino Mediterranean Observatory) project |161| is pursued by Italian 
groups and started 1998 with the objective to investigate the feasibihty of a cubic-kilometre- 
scale deep-sea detector and to identify and explore a suitable site. 

Technical solutions have been devised, investigated and optimised for an easy-to-deploy, 
stable and cost-effective detector |162H164] . The supposedly most important new concept 
developed in NEMO is that of "flexible towers" constructed of horizontal bars of a length 
of up to 15 m, interconnected by ropes forming a tetrahedral structure, so that adjacent 



bars are positioned orthogonal to each other (see Fig. 23). The optical modules are fixed 
to the bars. This design has several advantages, in particular (i) that a tower can be folded 
together and deployed to the sea floor as a compact object that is subsequently unfurled; 
(ii) the separation of the mechanical tension (carried by the ropes) from the electro-optical 
backbone cable; (iii) the provision of a 3-dimensional arrangement of photomultipliers per 
tower, allowing for local reconstruction of muon directions. 




Figure 23. Sketch of the prmciple of a flexible tower during the unfurhng process. 
Note that exact configuration and the packaging of the bars does not correspond to 
the current design framework. 



Concurrently with the technical developments, a suitable site at a depth of 3.5 km, about 
100 km off Capo Passero at the South-Eastern coast of Sicily has been identified and inves- 
tigated during various campaigns (see [103J and references therein). 

During the first prototyping phase ("NEMO Phase-1"), a cable to a test site near Catania 
at a depth of 2 km was installed and equipped with a junction box. In 2007, a "mini- 
tower" with 4 bars was deployed, connected and operated for several weeks. Although the 
data taking period was limited to a few months due to technical problems, the mini-tower 
provided the proof of concept for the technologies and most of the components employed. 
The flux of atmospheric muons was determined in good agreement with the expectations 
from simulation [165J. 

The Phase-2 [166J setup includes shore infrastructure at Capo Passero and a 100 km long 
cable to the site at 3.5 km depth; both are meanwhile in place. A remotely operated vehicle 
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(ROV) is available for the deep-sea operations. A mechanical test tower of limited size was 
successfully deployed and unfurled in early 2010. The plans to deploy a full-size prototype 
tower will be pursued in the KMSNeT framework. 

During the Phase-1 activities, an acoustic setup, Oz/DE (Ocean noise Detection Ex- 
periment), was operated for almost two years and yielded important data on the acoustic 
deep-sea background |167] and also on sperm whales |168| . 

4.6. NESTOR 

NESTOR [169] started out as a Greece-centred, originally international collaboration 
around the beginning of the 1990s. A series of workshops was held in 1991-1993 at the 
homonymous institute in Pylos at the West coast of the Peloponnesus and a detector concept 
presented in 1994 |170] . The basic idea is to attach the optical modules to rigid hexagonal 
structures ("floors") of 30 m diameter, of which 12 are stacked to a tower with 20-30 m inter- 
floor distance. The deployment would be from the sea top, the connection to the main 
cable or already-deployed components being made in the sequence recovery-connection- 
redeployment, thus avoiding the use of deep-sea submersibles. 

Several suitable site options off the coast of Pylos have been identified. They are between 
3000 m and 5200 m deep and have distances to shore of 30-50 km. The NESTOR collabora- 
tion has performed numerous studies on the site characteristics, most recently on the water 
optical properties (see |103| and references therein). 





Figure 24. Photograph of the Delta Berenike 
platform in the final phase of its construction. 

After a long phase of technical development a cable was installed to a site at 4 km depth. 
A single, reduced-size floor with 12 optical modules was deployed, connected and operated for 
more than a month; its operation had to be terminated due to a failure of the cable to shore. 
The data recorded sufficed to validate the detector functionality within specifications |171] 
and to measure the atmospheric muon fiux, which was found to agree with the expectations 

To facilitate the deployment of the fioors, NESTOR has developed and constructed 



a triangular working platform for sea operations, the Delta Berenike (see Fig. 24). The 
platform with side lengths of about 50 m is carried by three cylindrical pontoons at its 
apices and has a triangular aperture for deployment operations. 
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5. Second-Generation Neutrino Telescope Projects 



There are three projects which reach, or even exceed, the size originally conceived by 
the DUMAND pioneers: IceCube at the South Pole, in its basic configuration, has been 
completed in December 2010; KMSNeT (km^ Neutrino Telescope) in the Mediterranean Sea 
and GVD (Gigaton Volume Detector) in Lake Baikal are in their preparatory/prototyping 
phase. IceCube instruments Ikm^ of ice. GVD is also planned to cover Ikm^, but with a 
much higher energy threshold than IceCube. KMSNeT envisages an instrumented volume 
of several km^. KMSNeT and GVD could be completed by 2017 but will start data taking 
earlier with intermediate configurations. It is expected that these detectors will form a 
global network or even join to a global neutrino observatory. 

5.1. IceCube 

IceCube |17Sl I174j is the successor of AMANDA. It consists of 5160 digital optical mod- 
ules (DOMs) installed on 86 strings at depths of 1450 to 2450 m in the Antarctic ice |175| . 
S20 further DOMs are installed in Ice Top, an array of detector stations on the ice surface 



directly above the strings (see Fig. 25). AMANDA, initially running as a low-energy sub- 
detector of IceCube, was decommissioned in 2008 and replaced by DeepCore, a high-density 
sub-array of six strings at large depths (i.e. in the best ice layer) at the centre of IceCube. 
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Figure 25. Schematic view of the 
IceCube neutrino observatory. Also 
shown is the location of AMANDA (cf. 



Sect. 4.3) and of DeepCore, a nested 



low-threshold array (cf. Sect. 5.1.1). At 
the surface, the air shower array IceTop 



(cf. Sect. 5.1.2) and the IceCube counting 



house are indicated. The Eiffel tower is 
shown to scale for a size comparison. 



For IceCube construction, the thermal power of the hot-water drill factory has been 
upgraded to 5 MW, compared to 2 MW for AMANDA, thus reducing the average time to 
drill a 2450 m deep hole with a diameter of 60 cm to S5 hours. The subsequent installation 
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of a string required typically 12 hours. Installation of IceCube started in January 2005 with 
the first string and was completed with the deployment of the 86**^ string at Dec. 18, 2010. 

As the components are not accessible after refreezing of the holes, the IceCube architec- 
ture avoids single point failures in the ice. A string carries 60 DOMs, with 30 twisted copper 
pair cables providing power and communication. Neighbouring DOMs share the same wire 
pair and are thus connected to enable fast local coincidence triggering in the ice. 



A schematic view of a DOM is shown in Fig. 26 A 10-inch photomultiplier (Hamamatsu 
R7081-02) is embedded in a 13-inch glass sphere. A mu-metal grid reduces the influence of 
the Earth magnetic field. The programmable high voltage is generated inside the DOM. The 
average photomultiplier gain is set to 10^. Signals are sent to the main board where they 
are digitised by a fast analogue transient waveform recorder (ATWD, 3.3 ns sampling) and 
by a FADC (25ns sampling). The photomultiplier signal is amplified by 3 different gains to 
extend the dynamic range of the ATWD to 16 bits. The resulting linear dynamic range is 
400 photo-electrons in 15 ns; the dynamic range integrated over 2 |as is about 5000 photo- 
electrons [176j . The digital section of the main board is based on a field-programmable gate 
array (FPGA), which communicates with the surface electronics and allows for uploading 
new program code. LEDs on a "flasher board" emit calibration pulses at a wavelength of 
405 nm which can be adjusted in intensity over a wide range up to 10^^ photons. 
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Figure 26. Schematic view of an 
IceCube Digital Optical Module 
(DOM). The analogue 
photomultiplier signals are 
digitised on the main board. For 
details see |176) . 



All digitised photomultiplier pulses are sent to the surface. In order to compress data 
for isolated hits which are mostly noise pulses, the full waveform is only sent for pulses 
appearing in local (neighbour or next-to-neighbour) coincidences on a string. All DOMs 
have precise quartz oscillators providing local clock signals, which are synchronised every 
few seconds to a central GPS clock. The time resolution is about 2 ns. The noise rate for 
DOMs in the deep ice is about 540 Hz and is reduced to ca. 280 Hz if an artificial deadtime of 



50 1^ is applied to discard after-pulses (only for the supernova burst trigger, see Sect. 5.1.3) 



These very low noise rates are essential for the detection of the low-energy neutrino emission 



associated with a supernova collapse (see Sect. 5.1.3). 

At the surface, 8 custom PCI cards per string provide power, communication and time 
calibration. Subsequent processors sort and buffer hits until the array trigger and event 
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builder process is completed. This architecture allows for a deadtime-free operation. The 
raw data rate of the full array is about 800 GB / day, which are written to tape and processed 
online on a computer farm to extract interesting event classes, like up-going muon candidates, 
high-energy events, IceTop/IceCube coincidences, cascade events, events from the direction 
of the Moon or events in coincidence with Gamma Ray Bursts (GRBs). The filtered data 
stream (about 70GB/day) is transmitted to the Northern hemisphere via satellite. 

The muon angular resolution achieved with present reconstruction algorithms is about 
1° for 1 TeV muons. It is expected that improved reconstruction algorithms using the full 
waveform will lead to a resolution below 0.5° for energies above 10 TeV. The very clear ice 
below a depth of 2100 m has a particular potential for improved resolution. This will be even 
more important for the angular reconstruction of cascades, for which the presently achieved 
angular resolution is only 30° - much worse than for water, mainly due to the strong light 
scattering in ice. 



5.1.1. Deep Core 

The geometry of DeepCore is sketched in Fig. |27| DeepCore consists of 7 central stan- 
dard IceCube strings plus additional 6 special strings. Ten of the DOMs of these strings 
are arranged at 1750-1850 m depth, above the dust layer with bad optical transparency. 
They are used as veto-detector for the deeper component. The deep component comprises 
50 DOMs per string and is installed in the very clear ice at depths between 2100 and 2450 m, 
where the effective scattering length ranges up to 50 m and the absorption length to 230 m. 
The six additional DeepCore strings are equipped with photomultipliers with an enhanced 
quantum efficiency (Hamamatsu R7081-MOD). 
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Figure 27. Layout of the DeepCore 
sub-detector. Shown are the positions of the 
DeepCore DOMs; for a 3-dimensional 



representation see Fig. 25 The depth profile of 



the ice transparency is indicated on the left. 



DeepCore has a factor of about six better sensitivity in photon collection than IceCube, 
due to the smaller spacing between strings (72 instead of 125 m) and OMs along a string (7 
instead of 17 m); the better ice quality; and the 30% higher quantum efficiency of the new 
photomultipliers. Together with the veto provided by IceCube, this results in an expected 
threshold of less than 10 GeV. This opens a new window for oscillation physics not tested 
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by other experiments and allows for probing dark matter models not covered by direct 
searches. The veto will also allow for identification of neutrinos from above if they interact 
within DeepCore. This enlarges the field of view of IceCube to the full sky |177] . 

5.1.2. IceTop 

IceCube is the only large neutrino telescope which can be permanently operated to- 
gether with a surface air shower array, IceTop |178| . IceTop consists of tanks filled with 
ice, each instrumented with 2 DOMs. With IceTop, the energy spectrum of air showers can 
be measured up to primary particle energies around 10^^ eV. The combination of IceTop 
information (reflecting dominantly the electron component of the air shower) and IceCube 
information (muons from the hadronic component) restricts the mass range of the primary 
particle. The comparison of air shower directions measured with IceTop and directions of 
muons from these showers in IceCube establishes a tool for angular calibration of IceCube 
(absolute pointing and angular resolution). 

5.1.3. IceCube as a supernova burst detector 

Last but not least, IceCube can be operated in a mode that is only possible in ice: The 
detection of burst neutrinos from supernovae. The low dark-count rate of the photomulti- 
pliers (280 Hz, see above) allows for detection of the feeble increase of the summed count 
rates of all photomultipliers during several seconds, which would be produced by millions 
of interactions of few-MeV neutrinos from a supernova burst |179l 1180) . IceCube records 
the counting rate of all photomultipliers in millisecond steps. A supernova in the centre of 
the Galaxy would be detected with extremely high significance and the onset of the pulse 
could be measured in unprecedented detail. Even a SN1987A-type supernova in the Large 
Magellanic Cloud would provide a recognisable signal and be sufficient to provide a trigger 
to the Super Nova Early Warning System, SNEWS |181] . 

5.2. KMSNeT 

In recognition of the fact that at least a cubic-kilometre sized detector will be necessary 
to really observe abundant astrophysical high-energy neutrino sources, the High Energy Neu- 
trino Astronomy Panel (HENAP) of the PaNAGIC^ Committee of lUPAP^ has concluded 
in its 2002 report |182] that "a km^-scale detector in the Northern hemisphere should be 
built to complement the IceCube detector being constructed at the South Pole". 

Following this recommendation, the Mediterranean neutrino telescopes groups - together 
with deep-sea technology and marine science groups - have formed the KMSNeT collabo- 
ration to prepare, construct and operate such a device. In 2006-2009, a Design Study, 
supported with 9 M€ by the EU, was conducted; its major achievements are a Conceptual 
Design Report (CDR) [l83] and a Technical Design Report (TDR) [103]. A further EU 
project, the KM3NeT Preparatory Phase (2008-2012) provides resources and a framework 
for work directed towards solving the funding, governance, legal and strategic questions and 
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also for engineering activities. KM3NeT will be a deep-sea research infrastructure hosting a 
neutrino telescope, but also providing continuous, long-tern access to deep-sea measurements 
to a variety of science communities, such as marine biologists, oceanographers, geophysicists 
and environmental scientists. In recognition of its high scientific potential, KM3NeT has 
been included in the priority project list of the European Strategy Forum on Research In- 
frastructures, ESFRI pillTg^ . 

Based on the experience and expertise of the first-generation projects, a variety of new, 
cost-effective design solutions for the neutrino telescope have been elaborated. The original 
goal of reducing the capital investment for a cubic kilometre of instrumented sea water to 
200 M€ has been outperformed by a factor of at least 4. Whereas the TDR still presented 
a set of design options, in particular for the mechanical structure of the strings/towers 
("detection units (DUs)") and the optical modules, further investigations have meanwhile 
led to convergence on a specific design, which will be described in the following. Since final 
prototyping and deployment tests are still to come, the discarded solutions partly serve as 
backup options. 

For the mechanical DU structure, a solution along the NEMO design (see Sect. 4.5) has 
been chosen. The bars will be 6 m long and be equipped with one optical module (see below) 



at each end; a schematic view is given in Fig. [28| Oil-filled hoses in equipressure with the 
deep-sea water will be used as vertical cables for electrical leads and optical fibres. Two such 
backbone cables, wound loosely around the ropes, will be used per DU, with break-outs for 
two copper leads and one optical fibre at each optical module. The DUs will be stacked in 
a cubicle pile for transport and deployment, fitting into a standard shipping container. The 



deployment will proceed as discussed in Sect. 4.5 



Figure 28. Left: Schematic drawing of a bar of the KM3NeT DU. The overall length is slightly 
below 6 m. The red object in the centre is a buoy made of syntactic foam, the yellow cylinders house 
the ropes and vertical backbone cables before unfurling. Right: A full DU stacked for transport 
and deployment (top view). The round bluish object is the top buoy. 



The digital optical module (DOM, Fig. 29) will be a 17-inch glass sphere, equipped with 



31 3-inch photomultipliers, their high-voltage bases and the digitisation electronics. High- 
voltage bases with a power consumption as low as 140 mW for a complete optical module 
have been designed for this application. The photomultipliers are oriented from vertically 
downwards to about 45° upwards. They are supported by a foam structure and fixed to the 
glass sphere by optical gel. An aluminium structure is used to conduct the heat to the glass 
and to provide support for the electronic boards. Major advantages of this design are: 
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Figure 29. The KM3NeT digital optical module (DOM). Left: Photograph of a prototype; right: 
technical drawing, with the mushroom-shape heat conductor and the electronics components. 

• The overall photocathode area exceeds that of a 10-inch photomultiplier by more 
than a factor of three; a further increase is possible by extending the light collection 
area using reflective rings |186j . The number of penetrators and separate electronics 
containers, which are expensive and failure-prone, is thus reduced to a minimum. 

• Since the photomultipliers are read out individually, a very good one-vs.-two photo- 
electron resolution is obtained, which is essential for efficient online data filtering. 

• Some directional sensitivity is provided. 



The KM3NeT data acquisition will follow the all-data-to-shore concept (see Sect. 4.4). 
For each photomultiplier, the time intervals during which the analogue output signal ex- 
ceeds an adjustable threshold will be digitised and sent to shore. This time-over-threshold 
information allows for precise photo-electron counting for small signals and provides a loga- 
rithmic measure of the amplitude for large signals, i.e. effectively an infinite dynamic range. 
The technical implementation can either be achieved using a custom-designed ASIC [103J 
or a fast time-to-digital converter (TDC) on a field programmable gate array (FPGA). The 
latter option is new and might significantly simplify the system. The digitised data are 
transported to shore via optical point-to-point connections using the Gigabit ethernet pro- 
tocol and dense wavelength division multiplexing. The system is driven by on-shore lasers, 
which allows for easy maintenance and also for sending clock signals to the DOMs by the 
same laser beam that carries the photomultiplier signals on its way back to shore. 

Time and position calibration will be based on the same principles as in ANTARES (see 



Sect. 4.4). For the acoustic receivers on the DUs, the stand-alone hydrophones (that require 
extra penetrators) will presumably be replaced by piezo elements glued to the inner glass 
surface of the DOMs. 

The full KMSNeT neutrino telescope will consist of about 300 DUs. Their geometrical 
layout on the sea-fioor (the "footprint") is still subject to optimisation. For the sensitivity 
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studies in the TDR two equal-size homogeneous blocks of hexagonal layout were assumed, 
with an inter-DU distance of 180 m and a vertical distance between adjacent bars of 40 m. 
The partition into two blocks takes into account the facts that deployment and maintenance 
of the deep-sea cable network becomes increasingly difficult for large homogeneous setups, 
and that at least two main cables to shore will be needed for the full detector. 

Since the exact footprint is not yet known, the deep-sea cable network topology is still 
open. It will either be star-like with a primary and a set of secondary junction boxes 
or consist of a cable ring surrounding the detector, with several primary junction boxes 
connected directly to the DUs. In the cable network, single points of failure for the full or 
at least substantial parts of the detector are unavoidable, making recovery and maintenance 
operations necessary. In contrast, no maintenance is planned for the DUs themselves. 

The time-line of the further steps towards construction and operation of KMSNeT is 



shown in Fig. 30 The three sites of ANTARES, NEMO and NESTOR are under consider- 
ation for KMSNeT. The site question has scientific aspects that have been investigated in 
detail (depth, water transparency, bioluminescence etc., see [il03j and references therein); ad- 
ditionally, the site choice is strongly linked to the availability of regional funding sources and 
is therefore closely interwoven with a political decision on KMSNeT. It is currently hoped 
that a site decision can be made by late 2011 or early 2012. Subsequently, the final technical 
design will be worked out and presented in a detailed proposal. Assuming that funding, legal 
and administrative issues are sorted out by then, it will be possible to launch production 
at that point. Data taking will start as soon as the first DUs are operational. From a very 
early stage of its construction on, the data from the KMSNeT neutrino telescope will exceed 
data from first-generation Northern-hemisphere neutrino telescopes in quality and statistics 
and thus provide an exciting discovery potential. 



Design Study 
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Data taking 



CDR 
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Figure 30. Time-hne 
towards KMSNeT 
construction and operation. 



5.3. GVD in Lake Baikal 

The Baikal Collaboration plans the stepwise installation of a kilometre-scale array in 
Lake Baikal, the Gigaton Volume Detector, GVD |187] . It will consist of strings which are 



grouped in clusters of eight (see Fig. SI). This results in a relatively flexible structure. 



which allows for rearranging the clusters and meets best the deployment conditions from 
the ice. Each string carries 24 optical modules spaced uniformly from 900 m down to about 
1250 m depth |188| . The modules will house 10-inch photomultipliers (likely Hamamatsu 
R7081-HQE) with a peak quantum efficiency of about S5%. 
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Figure 31. The Baikal 
Gigaton Volume Detector 
(GVD). Left: 
Arrangement of the 12 
clusters; middle: One 
cluster; right: Schematic 
view of a string section. 



The optical modules on a string are grouped into two sections, each consisting of 12 OMs, 



a service module (SM) and the electronics unit (BEG) with its 200 MHz FADCs (see Fig. 31 



right). Analogue signals from the optical modules are transmitted to the BEG through 
coaxial cables. A trigger is formed by a coincidence of any neighbouring optical modules. 
Digitised waveforms for each triggered channel are transmitted via Ethernet from the BEG 
to the central underwater micro-PC of the cluster. The cluster DAQ provides inter-section 
time synchronisation, on-line data selection, and communication to shore through an optical 
cable. Prototype strings have been operated in 2009 and 2010. They have demonstrated a 
time accuracy of about 2 ns. In April 2011, a prototype cluster with three mini-strings and 
all key elements of DAQ electronics and the communication system was deployed. 

Simulations have been performed for 96 strings in 12 clusters and a total of 2304 OMs. A 
compromise between large volume for cascade detection and reasonable efficiency for muons 
was found for an instrumented height of 345 m, a cluster diameter of 120 m and a vertical 
spacing between optical modules of about 15 m [188]. At trigger level, the effective detection 
area for muons with energies above 3TeV is 0.2- 0.5 km^; the effective detection volume for 
cascades above 50TeV is 0.3- 0.8 km^. We note that cuts for background suppression will 
reduce these values significantly, in particular at lower energies. The directional accuracies 
are 0.5-1.0° for muons and 3-7° for cascades. 

A threshold for muons of about 3-lOTeV appears to be rather high when compared 
to IceCube and KM3NeT. On the other hand, the optimum energy cut to obtain the 
best signal-to-noise ratio (extraterrestrial versus atmospheric neutrinos) for the weakest 
detectable sources is at a few TeV for point sources, and in the 100 TeV range for diffuse 
fluxes |189] . Therefore, a sparse detector conflguration such as GVD may offer a favourable 
physics/cost ratio for neutrino signals extending to the 100 TeV region and beyond, as e.g. 
expected for Gamma Ray Bursts and AGN jets. 
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6. Physics Results and Perspectives 



After two decades of data taking with neutrino telescopes, a lot has been learned - 
despite the fact that not a single high-energy neutrino of cosmic origin has been clearly 
identified as yet. In this section, the results as of early 2011 are presented for atmospheric 



neutrinos (Sect. 6.1 ), for cosmic neutrinos including those from supernova bursts (Sect. 6.2) 



for searches for Dark Matter and other exotic particles (Sect. |6.3) and from cosmic-ray 



studies with IceCube/IceTop (Sect. |6.4) 



6.1. Atmospheric neutrinos 

Atmospheric neutrinos and muons are produced in cosmic-ray interactions in the atmo- 
sphere. Up to energies of about 100 TeV, their flux is dominated by pion and kaon decays. 
The corresponding neutrinos are referred to as "conventional" atmospheric neutrinos. The 
spectrum follows approximately an E~^''^ shape. It can be calculated with uncertainties of 
about 25% in the energy range lOOGeV-lTeV (see |190| and references therein, a widely 
used analytical calculation is given in |59]). At higher energies, "prompt" atmospheric neu- 
trinos from the decay of charm and bottom particles take over. These particles decay before 
having a chance for further interactions, and the resulting neutrinos therefore closely follow 
the primary cosmic ray power law spectrum, i.e. an shape. 

An almost background-free separation of neutrino-induced upward-going muons from the 
huge background of downward-going muons is the central requirement for an underwater or 



under-ice telescope (see Sect. 3.4.1). The first-generation experiments (Baikal, ANTARES 



and AMANDA) have quickly mastered this challenge, even more so IceCube. Figure 32 
shows the rate of muons as a function of the zenith angle 9 as measured with ANTARES. 
Below the horizon {9 < 0) the rate is well described by the expectation for atmospheric 
neutrinos, above the horizon by that for atmospheric muons. 




Figure 32. Number of reconstructed 
muons in the 2008 ANTARES data, as a 
function of the reconstructed zenith angle 
^Reco (black error bars). Also indicated 
are the simulation results for atmospheric 
muons (red dashed), and muons induced 
by atmospheric neutrinos (blue). The 
shaded band indicates the systematic 



'Reco uncertainties. Figure taken from |125) . 



The energy spectrum of atmospheric neutrinos is inferred from muon neutrino charged 
current events since they offer the best event statistics and, for this measurement, the least 
background. Measuring the energy spectrum is difficult since the weak dependence of the 
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Cherenkov light yield (i.e. of dE^/dx) on needs to be exploited; additionally, one has 
to take into account that most muons lose an unknown fraction of their energy before 
reaching the detector. The neutrino energy spectrum is therefore determined with the help 
of sophisticated deconvolution procedures, inducing significant point-to-point correlations. 
Two underground experiments have published atmospheric neutrino energy spectra up to 
a few TeV: The Frejus experiment |191] and Super-Kamiokande |192) . AMANDA and 
IceCube have extended this energy range by two orders of magnitude, up to 200 and 400 TeV, 



respectively. Figure [33] shows the spectra as published by these experiments. The data are 
well compatible with the predictions for conventional atmospheric neutrinos. In particular, 
no excess at high energies is observed as yet; improved data statistics from IceCube, however, 
will soon allow to test flux models for prompt neutrinos from the decay of charm and bottom 
hadrons, which would show up as a shoulder at some 100 TeV. In addition, constraints on 
the neutrino flux in this energy range and beyond are used to place upper limits on the flux 



of extra-terrestrial sources with a hard spectrum like E^^"^ (see Sect. 6.2). 
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Figure 33. Energy spectrum of 
atmospheric neutrinos. Green triangles: 
Frejus |191| : blue band: 
Super-Kamiokande [192j; red band: 
AMANDA forward folding analysis |98j; 
red triangles: AMANDA unfolding 
analysis |193| : black triangles: 
IceCube-40 unfolding analysis [194]. 



Atmospheric neutrinos also provide a tool to investigate neutrino oscillations. Standard 
oscillation lengths scale with E,^. For distances of the order of the Earth diameter the 
first oscillation minimum is at E^ ~ 24GeV (see Fig. 34). Violation of Lorentz Invariance 
(VLI), as suggested by certain quantum gravity theories, also leads to oscillation effects, 
with oscillation lengths that scale with l/E^, in the simplest case, n = I (see Sect. |2]). The 



left panel of Fig. 34 shows the survival probability as a function of neutrino energy for a 
baseline of the Earth diameter and assuming maximal mixing for all three cases (conventional 
oscillations, VLI oscillations and quantum decoherence, with the assumptions given in the 
figure caption). Large new neutrino telescopes are an ideal tool to look for non-standard 
oscillations, due to their ability to collect large statistics of high-energy neutrinos having 
travelled over distances of several thousand kilometres. From the non-observation of a 
deficit at high energies and the angular dependence, limits on the relevant parameters can 
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be derived. For n = 1 the 90% CL upper limits from AMANDA are 2.8 x lO"^? 
and 1.2 x 10^^'' for Di (see Sect. |2|. This is similar to limits obtained from the much 
smaller underground experiments. For n = 2 the limits are 2.7 X 10-31 GeV"^ for A5 
and 1.3 x lO^^i GeV"^ for D* = Di/E''-^ p]. IceCube promises another factor of ten in 
sensitivity over the forthcoming years. 




Figure 34. Left: survival probability as a function of energy for a baseline of the Earth diameter 
and for conventional oscillations (solid line), VLI (dotted line) with n = 1, A6 = 10"^^ and maximal 
mixing, and quantum decoherence (dashed line) with n = 2 and D* = 10"'^'^ GeV^^. See [98j for 
analysis methods taking into account the limited energy resolution. Right: allowed regions at 90%, 
95% and 99% confidence levels (from darkest to lightest) for VLI-induced oscillation effects (n = 1) 
from AMANDA data, together with the 90%-contour from a combined Super-Kamiokande and K2K 
analysis (dashed line) and the projected IceCube 10-year 90% sensitivity (dotted line). Figures are 
taken from |98| . 

Another exotic signature of VLI would be directional dependences of neutrino interac- 
tions or oscillations. Accelerator experiments like MINOS and K2K have searched for a 
sidereal modulation of their interactions rate (i.e. for a dependence on the orientation of 
their beam axes with respect to the Sun) and did not find any effect. A search for a sidereal 
effect of rates of atmospheric neutrinos in IceCube led to limits on the relevant parameters 
which improve MIN0S/K2K limits by factors of 3-1000, due to the longer baseline and the 
higher energy (see |195) and references therein). 

6.2. Cosmic neutrinos 

High-energy cosmic neutrinos may either be identified as accumulation of events pointing 
to a particular celestial direction ("point sources") or as extended diffuse emission, ranging 
from a few degrees (as for nearby supernova remnants) to fully diffuse, expectedly isotropic 
neutrino ffux; in both cases, the signal needs to be distinguishable from the intrinsic back- 
ground of atmospheric neutrinos. 
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Point-source searches use the directional and energy information to reduce this back- 
ground. Cosmic neutrinos from a given source would cluster around the source direction, 
with a point spread function determined by the angle between muon and neutrino and by 
the detector angular resolution for muons, both depending on energy. A further handle in 
point source searches comes from the fact that generic extraterrestrial sources have a harder 
spectrum than atmospheric neutrinos. For variable sources the time distribution may be 
used as additional criterion, in particular if independent information is available (as e.g. 
gamma observations in the case of GRB signals). 

Searches for diffuse fluxes can only use the measured energy as criterion for separating 
cosmic and atmospheric neutrinos. They thus critically depend on a detailed understanding 
of the detector response as a function of energy. Moreover, the high-energy tail of atmo- 
spheric neutrinos is dominated by prompt neutrinos, whose flux has larger uncertainties 
than that of conventional atmospheric neutrinos. 

In the following, we will first summarise the results on diffuse fluxes obtained so far, then 
those for steady point sources, and finally methods and results for variable sources. 

6.2.1. Searches for diffuse cosmic neutrino fluxes 

A diffuse neutrino signal may reveal itself as an excess of 

1. high energy upward moving muons from TeV-PeV muon neutrinos; 

2. high energy contained cascades; 

3. extremely high energy events emerging from downward moving PeV-EeV neutrinos 
of all flavours. 



The first case corresponds to the standard signature of muon neutrinos interacting via 
charged current reactions. Due to neutrino oscillations, the expected ratio of neutrino 



flavours at Earth is z/g : z/^ : z/^ = 1 : 1 : 1 (see Sect. 2.2.3). Since only charged cur- 
rent reactions of z/^ and, in about 17% of all cases, of z/^ produce high-energy final-state 
muons, whereas all other reaction channels appear as cascades, the relevance of detecting 
cascade events is obvious. The reduced directional accuracy (much worse for cascades than 
for muons) is not very important for diffuse fluxes. 

For first-generation detectors, the ratio between the "fiducial" volume (within which 
cascades can be reliably identified) and the effective volume for muon detection is small. 
Therefore the cascade channel does not necessarily provide superior sensitivity for diffuse 
fluxes. This is particularly true for detectors in ice where light scattering makes the iden- 
tification of cascades challenging. The situation is more favourable in water, where even 
events beyond the geometrical volume of the detector may be detected and reconstructed. 
Even without precise reconstruction of distant cascades, however, the non-observation of 
large signals from outside the geometrical volume can be used to derive upper limits on the 
diffuse flux. This approach has been pursued by the Baikal collaboration, who - in spite of 
the small volume of the NT200 detector - could compete with the much larger AMANDA 
detector over many years. 
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For detectors on the cubic-kilometre scale, with a much smaller surface-to-volume ratio, 
the situation is more favourable for the cascade channel. The relevance of contained cascades 
for diffuse flux measurements is enhanced by the fact that their energy deposit in the detector 
volume - and thus their Cherenkov light output - correlates much more strongly with the 
neutrino energy than that of muons. Furthermore, the outer part of the instrumented volume 
can be used as a veto against events mimicking cascades. 

No significant excesses over atmospheric neutrinos or other kinds of background has been 
observed so far, resulting in upper limits on the diffuse flux of extraterrestrial high energy 



neutrinos. Figure |35] summarises the limits obtained in the TeV - PeV region. For each 
experiment and each method only the best limit is shown. Remarkably, from the first limit 
derived from the underground experiment Frejus (1996) to the 2010 IceCube-40 limit, a 
factor of 500 improvement has been achieved. Several models as e.g. the blazar model of 
Stecker |57) shown in the figure can be excluded. A further factor of 10 improvement is 
expected over the next 2-3 years, using the full IceCube detector and combining muon and 
cascade information. The expected sensitivity is more than an order of magnitude below the 
Waxman-Bahcall bound, and prompt atmospheric neutrinos will be detectable for all but 
the lowest predictions |196] . Also, a test of the Waxman-Bahcall prediction for the diffuse 
flux of neutrinos from Gamma Ray Bursts [50j will be in reach. We note, however, that this 
model is much more easily tested with point source methods for transient, triggered sources 
(see below). 

Neutrinos in the multi-PeV to EeV region must come from above or from close to the 
horizon to be detected since the Earth is not transparent for neutrinos of such energies (see 



Fig. [9] in Sect. 3.1). The IceCube collaboration has performed analyses which are tailored 



to these neutrinos of extreme energies and yield limits |193[ 1203] that meanwhile are equal 
or better than those obtained from radio or air shower detectors (see Sect. jr)). Fi gure 36 



summarises the differential limits in the PeV -EeV range obtained by IceCube and by radio 
and air shower methods. 

The figure demonstrates that IceCube can merely detect cosmogenic neutrinos from 
cosmic ray interactions with the cosmic microwave background radiation. Also Auger (with 
a multi-year exposure) and ANITA (with another flight) will have little chances to detect 
this "guaranteed source". Only dedicated detectors like e.g. next-generation radio or acoustic 
detectors or maybe a very large, sparsely instrumented configuration of KM3NeT will provide 
the required sensitivities. Arguably, other sources than cosmogenic neutrinos may populate 
the highest-energy region, making explorative studies in this energy range valuable. The 
alternative detection methods (radio, acoustic, air shower detection) are addressed in Sect.[7j 

6.2.2. Searches for steady neutrino point sources 

In initial searches for point sources |209H211] . the sky was subdivided into bins of a size 
which was optimised to the detector resolution and the expected signal energy spectrum 
(typically E~'^). In these "binned searches", the signal would appear as an excess over 
atmospheric neutrinos in a certain bin. In order not to lose sources through signal sharing 
between adjacent bins, the search had to be repeated with shifted bins, resulting in trial 
factors which effectively reduced the sensitivity. Present searches, in contrast, use likelihood 
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Figure 35. 90% C.L. integral upper limits on the diffuse ffux of extraterrestrial neutrinos. The 
horizontal lines extend over the energy range which would cover 90% of the detected events from 
an E~'^ source (5% would be below and 5% above the range). All model predictions have been 
normalised to one flavour, i.e. all of the all-flavour limits have been divided by 3. The coloured 
band indicates the measured flux of atmospheric neutrinos (see also Fig. 33), the broadening at 



higher energies reflects the uncertainties for prompt neutrinos. The limits on muon neutrinos are 
from 807 days AMANDA [197|, 334 days ANTARES |198], and 375 days IceCube-40 (prehminary) . 
Cascade/all flavour hmits are from 807 days AMANDA [199], 1038 days Baikal-NT200 p^l200] . 
and 257 days IceCube-22 [201j. The Frejus and MACRO limits have been published in jl 1 j and 
|202j . respectively. Also indicated is the Waxman-Bahcall (WB) bound f60] . see Sect. 2.1.4 



functions which account for the smearing of the signal with a given point spread function. 
The probability of an event originating from a given source is calculated from a 2-dimensional 
probability density (typically a Gaussian), with the width defined by the uncertainty of the 
event direction. This "unbinned method" turned out to be up to 40% more sensitive than 
the binned |212) . Using an energy estimator to distinguish hard extraterrestrial from soft 
atmospheric neutrino spectra further enhances the sensitivity. 

The upper limit on a E~'^ neutrino flux from any point source in the Northern sky 
measured with AMANDA over seven years is £^^0 < 5.2 x 10~® GeV cm~^ s~^ sr~^ (averaged 
over the Northern hemisphere, using a binned method) |211) . Based on 813 days of data 
taking, ANTARES has recently released a sky map and preliminary values |213| for the 
flux limits of a set of preselected potential point sources (see Fig. 37), which improve the 
Southern sky limits from 14 years of Super-Kamiokande data. 

The real big leap, however, is being made with IceCube. First point source results 
were published for the 22-string configuration operated in 2007 |215| . In |216] this search 



53 




7 8 9 

log^„ (Energy/GeV) 

Figure 36. All-flavour 90% C.L. differential upper limits on the flux of extraterrestrial neutrinos 
in the PeV-EeV region. Limits are from the under-ice radio array RICE |204) . the air shower 
detector Auger [205] . the radio balloon experiment ANITA [206] . IceCube-22 [l93], and IceCube-40 
(preliminary) [203]. Also given are the expectations for 3 years of operation of the full IceCube 
detector and for a 1000 km^ shallow radio detector at the South Pole. The coloured band corre- 
sponds to different predictions for cosmogenic neutrinos from GZK processes. References for the 
GZK scenarios, the Z-burst prediction and the SUSY top-down scenario (TD) can be found in |203j . 
For recent predictions we refer to \64:\ I207| I208| . 



was extended up to 45° above the horizon. The background from down-going atmospheric 
muons was suppressed by hard cuts selecting only few thousand of the biUions of background 
events, including most of the PeV-EeV neutrinos. In that way, a reasonable sensitivity to 
the high-energy tail of hard neutrino spectra (e.g. was achieved. Naturally, the IceCube 
sensitivity to a corresponding source in the Southern hemisphere is worse than for Northern 
sources since the analysis relies exclusively on the tiny high-energy tail of the neutrino flux. 
Note, however, that IceCube would be almost blind to Southern sources with an energy 
cut-off in the PeV range or below (some remaining sensitivity is expected from DeepCore, 
which, using the rest of IceCube as a veto layer, could identify neutrinos which interact 
within the DeepCore volume with An acceptance). For unbroken E'"^ spectra, a cubic- 
kilometre detector at the South Pole can compete with a Northern first-generation detector 
like ANTARES up to a declination of 45°. This, on the other hand, means that there is a 
broad declination region where the combination of IceCube and ANTARES data will give a 
better sensitivity than IceCube or ANTARES alone. Such combined analyses are presently 
underway. 

Figure 38 shows the full-sky map derived from IceCube-40 data [214j taken between April 
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Antares 2007-2010, preliminary 




Figure 37. Left: Equatorial sltymap of neutrino-induced muon events from 813 days of ANTARES 
data from 2007-2010 (preliminary). The background colour scale indicates the sky visibility in 
percent of the time. The most significant accumulation of events, marked with a red circle, is 
compatible with the background expectation. Right: Flux limits for the 51 candidate sources 
marked in the skymap, compared to the expected average ANTARES sensitivity (blue dotted line) 
and to corresponding results from MACRO [209] . Super-Kamiokande | i210j . AMANDA [2TT] and 
IceCube [2Ti] . 



2008 and May 2009. The data sample contains 14 121 events from the Northern sky, mostly 
muons from atmospheric neutrinos, and 22 779 events from the Southern sky, mostly high 
energy atmospheric muons. With this sample, five searches were performed: (i) a scan of 
the entire sky for point sources; (ii) an analysis for 39 predefined potential sources (reducing 
the huge trial factor inherent to a full-sky search); (iii) the remaining three analyses stacked 
source candidates of different astrophysical nature and then searched for an excess. These 
classes were (iii-a) 16 sources of TeV gamma rays; (iii-b) 127 starburst galaxies and (iii-c) 
5 nearby galaxy clusters. The predefined sources included, among others, the Crab Nebula, 
the Geminga pulsar, the star cluster Cygnus 0B2 and the active galaxies Markarian 421, 
Markarian 501, 3C273 and M87. The results of the predefined point-source search (ii) are 



shown in Fig. 37 (right). All search results are consistent with the null hypothesis; the 



"hottest spot" from the all-sky scan is found at a right ascension of RA = 113.75° and a 
declination of DEC = 15.15° and has a chance value of 18%. Figure 39 compiles the limits 
from previous experiments, from the different IceCube stages and the expected sensitivity of 
one year of KM3NeT data. We note that the combined data of IceCube-40 and IceCube-59 
surpass the mark of 1 km^ x 1 year and thus exceed 1 year worth of data from the full IceCube 
detector. Very soon, a factor of 1000 improvement of the sensitivity to point sources will 
have been reached when compared to the very first AMANDA point source paper from 2000 

Where do we stand with respect to predictions for known sources? Figure 40 gives 
a "taste" of the answer. It shows the differential fluxes for three theoretical models and 
confronts them with the 90% C.L. upper limit and the 5cr discovery potential from IceCube- 
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Figure 38. Equatorial skymap 
of pre-trial significances of the 
all-sky point source search with 
IceCube-40 [214]. Each dot 
represents one neutrino event, 
the colour scale indicates the 
significance of event 
accumulations. The Galactic 
plane is shown as black curve. 
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Figure 39. Point source neutrino flux 
sensitivities (median expected limits at 
90% C.L.) from various experiments: 
Super-Kamiokande [210], AMANDA 
[2Tl] . IceCube-40 |2l5, IceCube-59, the 
sum of IceCube-40 and IceCube-59 
(preliminary results) and ANTARES 
|213| . Also shown is the expected 
sensitivity from one year of KM3NeT 
data |l03j . 



40. Its shows the z/^ + i?^ predictions for the supernova remnant RX-J1713. 7-3946 |i218j which 
was moved from its real position at the Southern sky the to the location of the Crab Nebula 
at the Northern sky, for MGRO J182+01 jSH], a gamma-ray source observed by the Milagro 
experiment, and for the active galaxy Centaurus A |219] . Note that recent measurements 
from the Fermi satellite [31j seem to disfavour a dominant hadronic origin of the gamma 
rays observed from RX-J1713. 7-3946 - the underlying assumption for the shown prediction 



see also Sect. 2.1.3). The conclusion is that optimistic model predictions are about one 



order of magnitude below present IceCube-40 limits. This does not rule out a discovery 
with a few years of full IceCube data. On the other hand, it seems that with IceCube - 
after having made the gigantic leap of a factor 1000 in sensitivity improvement - we are just 
scraping the discovery region. The main discovery potential therefore may remain for source 
phenomena which are not covered by the models addressed in Fig. 40, e.g. "dark" sources 



without significant high-energy gamma emission, or variable sources discussed in Sect. 6.2.3 



The implications for a Northern-hemisphere detector, however, are obvious: Its sensitivity 
must substantially exceed that of IceCube! 
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Figure 40. Differential ffux for 
tliree tlieoretical models of 
neutrino point sources, confronted 
with the 90% CL. upper limit and 
the 5a discovery potential from 
IceCube-40 (see text for 
explanations). Figure taken from 
[214j . 



6.2.3. Searches for variable neutrino sources 

Many astrophysical sources are known to have a variable flux at different wavelengths. 
Examples for such flaring sources are Active Galactic Nuclei, Soft Gamma Ray Repeaters, 
and Gamma Ray Bursts. Binary systems often show a periodic behaviour, as pulsars do. 
Neutrino searches for steady point sources integrate over time and continuously accumulate 
atmospheric neutrinos. That makes them less sensitive to the detection of flares when 
compared to searches over smaller time windows. Actually, compared to searches integrating 
over a full year, time-dependent searches can achieve 5a discoveries with two (flve) times 
less signal events for flares of 1 day (1 second) duration |220H222] - even if the there is no 
information on the flaring state from electromagnetic observations (so-called "untriggered" 
searches). In case this information is available and used ("triggered searches"), the gain can 
even be larger, in particular for Gamma Ray Bursts signalled by satellite detectors. 

IceCube has performed an untriggered search for flares from selected sources using 
IceCube-40 data (April 2008 to May 2009), and a triggered search for sources monitored by 
Fermi-LAT, SWIFT and Imaging Cherenkov Telescopes using IceCube-22 and IceCube-40 
data (May 2007 to May 2009) j223] . These analyses cover time intervals from 40 |js to a 
year. No signiflcant signal evidence was found. 

In this context it may be worth mentioning that of seven flares which were selected for 
the triggered search with IceCube-22, flve did not show a related neutrino signal in the 
selected time window, while one event was observed for each of the other two sources. One 
of these sources is the Active Galaxy lES 1959+650, which has a "neutrino history": In 
an analysis of AMANDA data taken from 2000-2003, flve events where recorded from the 
direction of 1ES1959+650. Interestingly, three of these came within 66 days in 2002. Two 
of the three neutrinos were coinciding within about a day with gamma-ray flares observed 
by the gamma-ray telescopes HEGRA and Whipple |224[ 1225] . Excitingly, one of these 2 



flares was not accompanied by an X-ray flare; such "orphan flares" would be expected for 
hadronic outbursts where the X-ray flux from synchrotron radiation of the electron plasma 
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is absent. This result was quickly followed by two theoretical papers, one claiming that 
the corresponding neutrino flux would not fit any reasonable assumption on the energetics 
of the source [226], the other claiming that scenarios yielding such fluxes were conceivable 
|227] . Since the analysis was not a fully blind analysis, it turned out to be impossible to 
determine chance probabilities for this event, and actually the result was never published in 
a journal. However, it initiated consideration to send alerts to gamma-ray telescopes in case 
time-clustered events from a certain direction would appear. These "Target-of-Opportunity" 
alert programs are described below. 

A special analysis [228J was performed for SN2008D, a core-collapse supernova at a 
distance of 27Mpc discovered by the SWIFT satellite on Jan. 9, 2008. Core-collapse su- 
pernovae might emit mildly relativistic jets (Lorentz factor F = 3-10) in which neutrinos 
could be produced by proton-proton collisions. Assuming that the jet pointed to the Earth, 
the non-observation of coinciding neutrinos with IceCube was used to constrain the total 
energy of the released energy and F. Actually, according to current models, the full IceCube 
detector could detect up to 100 events for a core-collapse supernova at lOMpc distance. 

One of the most promising sources of high energy neutrinos are Gamma Ray Bursts 



(GRBs). As mentioned in Sect. 2.1.4, neutrino emission can be modeled for three GRB 
phases: the precursor phase when the jet is still forming and no electromagnetic radiation is 
escaping [51]; the prompt phase coinciding with the burst in gamma rays (see e.g. [501 1229] ); 
and the afterglow phase [52j. Direction, time, duration and gamma spectrum of GRBs 
are provided by satellite observations. Early searches for neutrino events coinciding with 
GRBs have been performed with Super-Kamiokande, Baikal and AMANDA. The AMANDA 
analysis [230J looked for neutrinos coinciding with any of 408 well-located GRBs recorded 
in the Northern hemisphere between 1997 and 2003. No neutrinos were found during or 
immediately prior to the GRBs. An upper limit on the diffuse flux from all GRBs was 
derived, which was still a factor of 1.5 above the flux predicted in [50] for the prompt phase, 
but already two times below the precursor model published in [51] . In a recent analysis, 
neutrinos recorded with IceCube-40 have been analysed with respect to coincidences with 
any of 117 Northern- hemisphere GRBs recorded in the time interval April 2008 to May 2009 
[231] : preliminary results for IceCube-59 are also already available [232J. Differently to the 
AMANDA analysis, the expected neutrino spectra have been calculated according the the 
observed gamma spectra from each individual GRB and following the prompt phase model 
of [229) . Again no excess of neutrinos close in time and direction is observed, resulting in 



a limit below the predictions of [50| 1229] . see Fig. 41 Even though this does not rule out 
the general picture of the GRB fireball model of [50j which assumes that cosmic rays of 
highest energy essentially all emerge from GRB, the amount of energy of the accelerated 
protons transfered to pions is obviously smaller than assumed in these models. The transfer 
parameter is uncertain by a factor of 3 or more; this and other uncertainties in the model 
parameters leave room for a possible confirmation of these models with future IceCube data. 

6.2.4- Alert programs 

Neutrino telescopes monitor essentially a full hemisphere. In contrast, most gamma-ray. 
X-ray and optical observatories can observe only a small fraction of the sky at any given 
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Figure 41. 90%C.L. upper limits on the neutrino fluxes from GRBs set by AMANDA, IceCube- 
22, IceCube-40 and IceCube-59 (preliminary). The IceCube limits are calculated with respect to 
the flux expected from the model of Guetta et al. |229) . the AMANDA limit with respect to the 
Waxman-Bahcall flux [SU] which assumes an average shape of the GRB spectra. The fluence limits 
for single GRBs have been translated to a limit on the diffuse neutrino flux from GRBs, multiplied 
with E"^. The corresponding model predictions are indicated as dashed lines. The strikingly worse 
limit of IceCube-22 as compared to AMANDA and later IceCube configurations is due to small 
statistics and one observed conicidence. Figure adapted from [232]. 



moment. Therefore only a small subset of neutrino data can be used for correlation studies 
with these data. The ability to identify such correlations can be improved by optical follow- 
up (OFU) programs |233H235] . In these programs, the observation of doublets or triplets 
of neutrino events from the same direction (or one exceptionally energetic single neutrino 
event) triggers a small network of automated 1-2 metre telescopes to point in that direction. 
If the neutrinos are emitted from jets in core-collapse supernovae in other galaxies or from 
GRBs, the optical telescopes could identify the rising light emission from the supernova or 
the GRB afterglow. As was shown in [233J , OFUs can improve the sensitivity to neutrinos 
from supernovae and GRBs by a factor of 2-3. Both IceCube and Antares are running 
OFU programs, triggering the optical telescopes TAROT and ROTSE (for details see |235| . 
No coincidences have been reported so far, although the IceCube analysis provided some 
interesting (but not significant) doublets. Presently, these programs are extended to include 
the X-ray satellite SWIFT. 

Another follow-up program in IceCube is known as Neutrino Target of Opportunity 
(NToO) |236) . It was motivated by the coincidence in 2002 between two neutrinos from 



the direction of the AGN 1ES1959+650 and gamma ray flares from the same source (see 
Sect. 6.2.3). In this case, the problem is complicated by the comparatively long duration of 
AGN flares (hours to weeks instead of seconds to minutes for supernovae or GRBs). Doublets 
within such long time intervals appear too often from any direction of the sky to be useful for 
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such a program. Therefore only selected sources, which are known to show flaring behaviour, 
are used for the NToO. In case a "cluster in time" of events from the direction of one of these 
sources has accumulated, a trigger is sent to the MAGIC gamma-ray telescope in La Palma. 
As soon as possible MAGIC is then pointed to that direction. 

A further spectacular alert program is the "standard" IceCube supernova alert. As men- 



tioned in 5.1.3 , the low dark noise rate of the photomultipliers allows for a mode of operation 
which is best suited for ice: the detection of burst neutrinos via the feeble increase of the 
summed count rates of all photomultipliers. This increase would be produced by millions of 
neutrino interactions at energies up to some tens of MeV within several seconds |179l 1237] . 
The tiny amount of light from one of these interactions would usually fire only the closest 
photomultiplier, so that no event-by-event reconstruction is possible. With this method, a 
supernova in the centre of the Galaxy would be detected with extremely high confidence, and 
the onset of the pulse could be measured in unprecedented detail since IceCube records the 
counting rate in millisecond steps. Even a SN 1987A-type supernova in the Large Magellanic 
Cloud would be identifiable. 

Since neutrinos leave a supernova a few hours before light is emitted, neutrino signals 
can be used to issue an early alert to optical astronomers. Actually, IceCube alerts are fed 
into the SuperNova Early Warning System, SNEWS [181]. Currently, the detectors Super- 
Kamiokande (Japan), LVD and Borexino (Italy) as well as IceCube contribute to SNEWS, 
with a number of other neutrino detectors and gravitational wave detectors planning to join 
in the near future. 

The average supernova rate for a galaxy like ours is estimated to be 2-5 per century 
(see |238j for an overview). Since 1987 no supernovae have been observed in our Galaxy 
or its neighbourhood. There might have been, however. Galactic supernovae obscured by 
matter between them and us which consequently could have been visible only in neutrinos. 
The non-observation of low-energy neutrino signals by underground neutrino detectors like 
Super-Kamiokande and by AMANDA/IceCube sets a 90% C.L. upper limit of about eight 
per century on the average rate of such bursts (we have rescaled the 2007 limit obtained in 
|238] for 25 years of non-observation to meanwhile 29 years). 

6.3. Dark matter and other exotic particles 

In addition to neutrinos from astrophysical objects, as discussed above, neutrino tele- 
scopes are also sensitive to potential neutrino fluxes from dark matter annihilations (indirect 



dark matter searches) and to hypothesised exotic particles (see Sect. 2.2). The status of the 
corresponding experimental search results and the future expectations are discussed in this 
section. 

6.3.1. Dark matter 



As explained in Sect. 2.2, dark matter could accumulate in the Galactic centre, in the 
Sun or in the centre of the Earth and eventually release neutrinos via self-annihilation 
processes. These neutrinos would be detectable as an excess over the irreducible background 
of atmospheric neutrinos. Several underground and underwater /ice experiments (Baksan, 
MACRO, Super-Kamiokande, Baikal, AMANDA, IceCube (see |239] and references therein) 



60 



as well as ANTARES |24UH242| ) have searched for an excess of neutrinos from the centre 
of the Earth or the Sun, IceCube also for signals from the Galactic halo and the Galactic 
centre |243l 1244) - all without having identified any significant signal. Alternative indirect 
searches using gamma rays [245J or other particle messengers |i246j have not found clear dark 
matter signals either. 



0.05 <Slj h- <0.20 

^ ■ IC/AMANDA 2001-2006, bb 

• ■ IC/AMANDA 2001-2008, W" 



M5SM muon flux limits 



IceCube Preliminary 

tr5l<tr;i" CDMS<2010)+XENON100(2011) 
IC36 lao days sensitiuity, ir^W" '■ 
Super-K2011, bb 
Super-K2011, W~W- 




10' w 10' 

WIMP mass (GeV) 

0.Q.^<1^-^^' <Q.20 M55M 5D cross-section limits IceCube Preliminary 

trsi itrj"' CDMS(2010)-^XENON100(2011) Super-K 2011, ■ 

■ . IC/AMANDA 2001-2008, Super-K 2011, H'" ir _ 

■ ■ IC/AMANDA 2001-2008, H-»^ KIMS 2007 

, , COUPP 2008 

IC86 180 days sensitivity, W'"'nr 

Picasso 2009 

\ '■■ t*t" for mj < mj, CDMS 2010 " 




WIMP mass (GeV) 



Figure 42. 90% C.L. exclusion 
limits on the muon flux induced by 
neutrinos from neutralino 
annihilations in the Sun (top) and 
on the spin-dependent 
neutralino-proton cross section 
(bottom) from indirect dark matter 
searches in the neutrino channel, as 
functions of the WIMP mass. The 
shaded areas indicate the MSSM 
parameter regions not excluded by 
direct searches (see text). Indicated 
are the limits obtained with 
AMANDA, IceCube-22 and 
IceCube-40 (squares) assuming soft 
(66) and hard (VF+I^~) neutrino 
spectra from WIMP annihilations; 
the expected sensitivities of IceCube 
with DeepCore (red lines); and 
exclusion limits obtained by other 
experiments. Limits have been 
rescaled to a common muon 
threshold of 1 GeV. See [211] for 
references to the various 
experimental data; the 
Super-Kamiokande results are taken 
from [247]. 



The most sensitive search for neutrino-induced muons from the direction of the Sun has 
been performed with IceCube-22, using data taken between March and September 2007 
|244] ■ In these months, the Sun is below the horizon at the South Pole. Simulations have 
been run for several neutralino masses each with two annihilation channels, a hard chan- 
nel (annihilation into I^"*"!^"), and a soft channel (annihilation into bb). Figure 
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(left) 

shows the limits on the flux of muons induced by neutrinos produced in MSSM neutralino 
annihilations in the centre of the Sun as a function of the neutralino mass m^. Since the 
average energy of the muons would increase with neutralino mass and hardness of the de- 
cay channel, and since the muon detection probability increases with energy, the limits are 
tighter for harder decay spectra and for larger m^. The green shaded area represents flux 
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predictions from presently allowed combinations of the MSSM parameters based on direct 
search experiments. The grey shaded area represents the parameter space that would still 
be allowed if direct searches would improve their sensitivity by a factor of 1000. 

Since the muon flux is proportional to the capture rate of neutralinos in the Sun, the 
muon flux limit can be converted into a limit on the neutralino-proton scattering cross 
section, in particular its spin-dependent part, cr^+p (the Sun is mainly a proton system). 
In |244) . this conversion has been performed by assuming equilibrium between capture and 
annihilation rates. A conservative result is obtained by further assuming that the spin- 
independent cross section vanishes. The result of this analysis is given in Fig. 42 (right) and 
compared to other limits from indirect searches by Super-Kamiokande and IceCube-22 as 
well as the best current limits from direct searches. The figure illustrates the potential of 
indirect searches to explore the parameter space of spin-dependent neutralino interactions 
with a sensitivity exceeding by far that of direct searches. 

Using the same method, data from AMANDA and IceCube have been used to set limits 
on other dark matter candidates, like the lightest Kaluza-Klein particle (hypothesised in 
theories with extra dimensions, with an expected mass range similar to that for neutralinos) 
and Simpzillas (super-heavy relic particles with masses in the range of 10^-10^^ GeV). We 
refer to pH ES] for details. 

An alternative approach was pursued to estimate the KMSNeT sensitivity to neutrinos 
from dark matter annihilations. The mSUGRA parameter space was scanned; for each 
parameter set the resulting neutrino spectrum was calculated and the resulting event rate 
in KMSNeT determined, taking properly into account neutrino propagation through the 
Sun. The resulting limits for the spin-dependent cross section, equivalent to Fig. 42 (right), 
are shown in Fig. 43 (see [103J for further explanations and references to the experimental 
data). KMSNeT in the currently planned configuration (see Sect. 5.2), with its comparably 
high detection threshold, is less sensitive to indirect dark matter signals than IceCube with 
DeepCore. 
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Figure 43. Spin-dependent 
WIMP-proton cross section as a function 
of the neutralino mass. The coloured 
points correspond to mSUGRA 
parameters consistent with constraints 
from cosmic microwave measurements. 
The models within (beyond) the 
discovery reach of KM3NeT are coloured 
green (red). Limits of selected other 
experiments are also indicated. The 
KM3NeT and IceCube limits are for 10 
years of data taking. 
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6.3.2. Magnetic monopoles, Q-balls and nuclearites 

In water and ice detectors, relativistic magnetic monopoles can be identified via their 
strong light emission, 8300 times more intense than that of a minimally ionising muon. 
Even below the Cherenkov threshold (velocity Vc ~ 0.75c), down to v ^ 0.5c, the light 



emission is large due to accompanying delta electrons. Figure 44 shows the monopole flux 
limits obtained by Baikal, AMANDA and ANTARES, where only in the ANTARES anal- 
ysis light emission below fmon = 0.75 c was addressed. The Baikal flux limit of 4.7 x 
10"^'' cm~^ s~^ sr~^ for upward-moving monopoles with /3 ^ 1 is based on 1040 live days of 
NT200 (1998-2002) |2l9]; that of AMANDA of 2.8 x 10"^^ cm'^ s'^ sr'^ for upward-going 



and 3 X 10 cm ^ s ^ sr ^ for downward-going monopoles on the data taken with AMANDA 
in a single year (2000) |250| : the preliminary ANTARES limit for upward- moving monopoles 

-1 [251] 
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cm 
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the preliminary IceCube limit derived from 
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one year of data taking with the 22-string configuration is about 3 x 10~ cm 
|252] . IceCube in its full configuration will be yet another factor 3-4 more sensitive and thus 
advance to more than three orders of magnitude below the Parker bound [253j, which is set 
by the condition that magnetic monopoles must not destroy the Galactic magnetic fields. 
Note that due to the absence of predictions for magnetic monopole fluxes a non-observation 
cannot exclude certain theoretical models (as it is the case for dark matter searches). 
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Figure 44. 90% C.L. upper limits on the flux 
of relativistic magnetic monopoles as a 
function of their velocity, /^mon = "Wmon/c- The 
ANTARES and IceCube results are 
preliminary. 



Slow GUT monopoles, Q-balls and nuclearites (see Sect. 2.2) would heavily ionise the 



matter they traverse, or even catalyse proton decays along their path. Such particles with (3 
in the range 10~^ - 10~^ have been searched for with underground and underwater/ice detec- 
tors; in absence of any signals, upper limits on their fluxes of the order of 10~^^ cm~^ s~^ sr~^ 
are derived [.254--257J . Note that limits based on the catalysis of proton decay are conditional 
since they depend on an assumed catalysis cross section. For instance, the limit of the Baikal 
experiment on GUT monopoles is 2[50] x 10^^^ cm^^ s^^ sr^^ for catalysis cross sections of 
10~^^[10~^''] cm^ and /3 = 10~^. Again, IceCube data will improve this sensitivity by about 
two orders of magnitude. Needless to mention that intense light emission of a slow object 
would be a spectacular signature which would make the first clear observation of such an 
event a real sensation - regardless of its concrete interpretation as GUT monopole, Q-ball 
or nuclearite. 
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6.4- Cosmic ray physics 

In detectors underwater/ice and underground, cosmic rays are recorded via the punch- 
through muons generated in air showers above the detector. Ideally, the deep detector 
is complemented by a surface air shower array for such measurements. In this case, the 
air shower array essentially responds to the total energy contained in the electromagnetic 
shower component, whereas the deep detector responds to the hadronic shower component 
from which the muons originate. Combination of both pieces of information allows for dis- 
entangling energy spectra and mass composition of the primary particles on a statistical 
basis. This approach has been pioneered with underground/air shower hybrid detectors 
like Baksan/Carpet (Russia) or MACRO/EAS-Top (Italy), followed by AMANDA/SPASE 
(South Pole Air Shower Experiment) and recently IceCube/IceTop (both at South Pole). 
The experimental results indicate that the average mass of the primary particles increases 
with energy in the knee region around 10^^ eV, see |258] and references therein. The com- 
bination IceCube/IceTop covers an enormously enlarged solid angle and an area of a full 
square kilometre, and yields data of unprecedented quality. With its huge statistics it will 
extend the energy range up to 10^^ eV, where the transition from Galactic to extragalactic 
cosmic rays is expected. 

Surface detectors can also be used to estimate the angular resolution of the deep detec- 
tor, by comparing the reconstructed directions of air showers and the corresponding muon. 
Comparing data from SPASE and AMANDA, the muon angular resolution of AMANDA 
was measured to be 2° -2.5° |259| . The angular resolution of IceCube was determined via 
the shadow of the moon with respect to cosmic rays. This shadow can be measured via 
down going muons with a high significance within a few weeks, and confirms that angular 
resolution and absolute pointing of IceCube are about 1° or better |260] . 

One of the initial motivations for deep-water detectors was to clarify puzzles in the depth- 
intensity relations that had emerged from underground muon measurements (see Sect. 4.1). 
Due to the energy dependence of the muon range in matter, these measurements are sensitive 
to the energy spectrum of muons from cosmic-ray-induced air showers. Figure 45 [26 Ij shows 
a compilation of recent data, dominated by muons with energies below 100 TeV; the measured 
muon flux excellently follows a model prediction [262J, thus demonstrating that this aspect 
of cosmic ray physics is well under control in the energy range covered. 

Over the last years, however, a new cosmic ray puzzle has emerged: Several detectors 
located at the Northern hemisphere have found tiny anisotropies in the arrival direction of 
cosmic rays in the range of several TeV to several hundreds of TeV (Milagro |263] . Tibet 
air shower array |264) . ARGO/YBJ [265j ) and of down-going atmospheric muons (Super- 
Kamiokande |266] ). Using recorded data of several 10^° down-going muons with energies 
above a few TeV, IceCube has detected per-mille anisotropies on all angular scales down to 
about 15°. A multipole analysis yields the strongest power for dipole (/ = 1) and quadrupole 
(/ = 2) contributions. The dipole does not point to the direction suggested by the Compton- 
Getting effect (anisotropy induced by the movement of the Earth [267J) if the cosmic ray 
plasma were at rest with respect to the Galactic centre. This indicates that Galactic cosmic 
rays co-rotate with the local Galactic magnetic field. Also smaller structures on scales 
between 15° and 30° are visible. The relative amplitude of the smaller structures is about a 
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Figure 45. Depth-intensity relation 
for muons measured by various 
detectors underground and 
underwater /ice. The shaded band 
indicates the systematic uncertainty 
of the ANTARES analysis, the 
dashed curve is the expectation 
derived from |262) . Note that single 
experiments can cover extended 
depth ranges by investigating muons 
incident under different zenith 
angles. Figure taken from |26T] . 



factor of five weaker than that of the dipole/quadrupole structure. The minima and maxima 
on all scales are reproducible in all hitherto analysed data samples, from IceCube-22 
to IceCube-59 [269]. 



Figure 46 shows the combined Milagro/IceCube significance sky-map after applying 



a band-pass filter to remove small-scale fluctuations (10° smoothing) and the dipole and 
quadrupole terms. 



Milagro + IceCube TeV Cosmic Ray Data (lo" Smoothing) 




-12 



-4 4 

significance 



12 



Figure 46. Significance skymap 
of the medium-scale anisotropies 
of atmospheric muon arrival 
directions in IceCube, combined 
with the Milagro cosmic-ray 
sky-map for the Northern 
hemisphere. The colour code 
indicates the significance of 
deviations from the average in 
a. Note that structures at very 
small and very large scales have 
been filtered out (see text). 
Figure taken from |269) . 



The Milagro data contain 2.2 x 10^^ air showers with a median energy of about 1 TeV, 
the overall IceCube data contain 3.4 x 10^° muons with a median energy of the primary 
particle of about 20 TeV. The origin of the medium scale anisotropies is still unknown. 
Possible explanations include a magnetic "nozzle" to a nearby supernova remnant, with a 
scatter-free transport along a field connected to the source |270| : at lower energies and small 
scales, the heliotail may play a role |271j . What ever the origin of the anisotropies is: These 
structures observed with IceCube may contain key information on nearby sources and/or 
cosmic ray propagation. To understand the anisotropies is an interesting challenge that will 
be re-addressed with the increasing statistics expected in the coming years. 
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7. Alternative Detection Principles for Extreme Energies 



The alternative detection technologies described in this section are tailored to signals 
which propagate with kilometre-scale attenuation. Consequently, they allow for the obser- 
vation of much larger volumes than those achievable for optical neutrino telescopes. Detec- 
tors on the 100 to 1000 km^ scale are necessary, for instance, to record more than just a 
few cosmogenic neutrinos, in the typical energy range of 100 PeV to lOEeV. This section 
is intended to give an outlook to future experimental opportunities. Some results of initial 
projects are included in Sect. [6} see in particular Fig. 36 



7.1. Detection via air showers 

At energies above 10^'' eV, large air shower arrays like the Pierre Auger Observatory 
(PAO) in Argentina [272] or the Telescope Array in Utah, USA [273j are seeking for hori- 
zontal air showers induced by neutrino interactions deep in the atmosphere (showers caused 



by charged cosmic ray interactions start much higher up in the atmosphere). Figure 47 
explains the detection principle. PAO consists of an array of water tanks covering an area 
of 3000 km^ that record the Cherenkov light of charged air-shower particles crossing the 
tanks. The array is combined with telescopes recording the atmospheric fluorescence light 
from air showers. The optimum sensitivity window for this method is at 1 - 100 EeV, the 
effective target mass is up to 20 Gigatons. An even better sensitivity might be obtained for 
tau neutrinos, z/^, scratching the Earth and interacting close to the array |274l 1275) . The 
charged r lepton produced in charged-current interaction can escape the rock around the 
array (in contrast to electrons) and mostly decays into hadrons (branching ratio ca. 65%) 
after a short path length (in contrast to muons). If this decay happens in the fleld of view of 
the fluorescence telescopes, the decay cascade can be recorded. Provided the experimental 
pattern allows for a clear identiflcation, the acceptance for this kind of signals can be large. 
For the optimal energy scale of 1 EeV, the present Vj. limit for an E'"^ tau neutrino flux is 
about E'^cf) < lO-^GeVcm'^g-ig^-i 

Space-based observation of extended air showers is an approach to even further increase 
the target mass for highest-energy cosmic-ray and neutrino detection, at energies beyond 
10^^ eV. The "Extreme Universe Space Observatory (EUSO)", proposed around the year 
2000, is a wide- fleld camera observing the atmosphere from an orbit at several 100 km height 
and registering the fluorescence light from extended air showers and, if circumstances allow, 
also the reflection at Earth surface of the Cherenkov light emitted in shower direction |277] . 
After several years of technical development and a long phase of uncertainty concerning the 
space carrier, plans are now to install the device - meanwhile renamed to JEM-EUSO |278] 
- on the Japanese Experiment Module (JEM) of the International Space Station, with the 
launch expected around 2015. JEM-EUSO will observe an atmospheric target volume with 
a mass of more than one Teraton and will thus exceed the PAO sensitivity by two orders 
of magnitude for energies above some 10^^ eV |279| . Similarly as for PAO, neutrino-induced 
air showers can be separated from those from hadrons or gammas by the depth of the 
interaction point in the atmosphere. The JEM-EUSO physics opportunities in the neutrino 
channel have been studied in general |280l 1281 ] and speciflcally for cosmogenic neutrinos 
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and GRB neutrinos 12831. It is not obvious that JEM-EUSO will detect neutrinos 



from known astrophysical or cosmogenic sources, but its measurements will explore possible 
neutrino fluxes in a hitherto inaccessible energy region. 

1.2. Radio detection 

Electromagnetic cascades, e.g. generated by high energy neutrino interactions in ice or 
salt, emit coherent Cherenkov radiation at radio frequencies. The effect was predicted in 
1962 |284] and confirmed by measurements at accelerators |285[ 1286) . Electrons from the 
material traversed are swept into the developing shower, which thus acquires an electric 
net charge. This charge propagates like a relativistic pancake of about 1 cm thickness and 
10 cm diameter. For wavelengths exceeding the cascade diameter, coherent emission of 
electromagnetic radiation (Cherenkov radiation and synchrotron radiation caused by the 
Earth magnetic field) occurs. The signal amplitude increases with the square of the net 
charge in the cascade, i.e. it is proportional to E^, thus making the method particularly 
attractive for high-energy cascades. The resulting bipolar pulse is in the radio frequency band 
and has a width of 1 - 2 ns. In ice, attenuation lengths of more than a kilometre are observed 
for radio signals, depending on the frequency band and the ice temperature, implying that 
for energies above a few 10 PeV radio detection becomes competitive or superior to optical 
detection (with its attenuation length of the order 100m) [287]. 

The prototype "Radio Ice Cerenkov Experiment (RICE)" was operated at the South 
Pole, with 20 receivers and emitters buried at depths between 120 and 300 m. From the 
non-observation of very large pulses, limits on the diffuse flux of neutrinos with E > 100 PeV 
and on the flux of relativistic magnetic monopoles have been derived |288 



The "Antarctic Impulsive Transient Array (ANITA)" [289j is an array of radio antennas 
which has been flown at a balloon on an Antarctic circumpolar path in 2006 and 2008/09 



see Fig. 48 , left). From 35 km altitude it searched for radio pulses from neutrino interactions 



in the thick ice cover and monitored, with a threshold in the EeV range, a volume of the 
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order of 10^ Gigatons. The resulting neutrino flux limits [206j are presented in Fig. 36 in 
Sect, ion 



A acoustic 




Figure 48. Principles of radio and acoustic detection of neutrinos. Left: ANITA balloon experi- 
ment. Right: Acoustic emission of a particle cascade. 



Future plans for radio detection in ice foresee large arrays of antennas on the surface of 
the antarctic ice shelf (ARIANNA [290J) or in the South Polar ice close to the surface (ARA 

m)- 

Even more exotic is the search for radio emission from extremely-high energy cascades 
induced by neutrinos or cosmic rays skimming the moon surface. An example is the "Gold- 
stone Ultra-high Energy Neutrino Experiment (GLUE)" which used two NASA antennas and 
reached a maximum sensitivity at several ZeV = 1000 EeV |292] . With the same method, 
the NuMoon experiment at the Westerbork Radio Telescope was searching for extremely en- 
ergetic neutrinos |293] . Similar activities are under preparation in the context of the LOFAR 
experiment 



7.3. Acoustic detection 

Production of pressure waves by charged particles depositing energy in liquids or solid 
media was predicted in 1957 |295] and experimentally proven with high-intensity proton 
beams two decades later j296] . In the case of a particle cascade, its entire energy is deposited 
into the medium, mostly through ionisation, and converted to heat on a time scale that is 
very short compared to the typical time scales relevant for generation and propagation of 
acoustic pulses. The effect is a fast expansion, generating a bipolar acoustic pulse with a 



width of a few ten microseconds in water or ice (see Fig. 48, right), corresponding to a peak 
signal power at 20 kHz. Transversely to the pencil-like cascade, the acoustic pulse propagates 
into the medium within a disk-shaped volume with a thickness corresponding to the cascade 
length of about 10 m. Exploiting this method would require to detect the acoustic pulses 
on the background of the ambient and intrinsic noise with a sparsely instrumented detector. 
This implies a very high detection threshold, in the EeV range at best. 
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Acoustic detection is an option both for ice and sea water. For ice, the signal itself is 
expected to be higher and ambient noise to be lower than in sea water. A test array, SPATS 
(South Pole Acoustic Test Setup), has been deployed at the South Pole in order to determine 
the depth dependence of the speed of sound |297) . the attenuation length of acoustic signals 
and the ambient noise [299J . The results of the latter two measurements are slightly 
discouraging as the attenuation length turns out to be about 300 m, an order of magnitude 
smaller than expected, and the noise level is roughly the same as in the deep sea at a calm 
sea state. As a variation of the ice approach, even the use of permafrost as medium has 
been discussed [300J. 

Test of acoustic detection in sea water have been performed close to Sicily (Oi^DE setup. 



see Sect. 4.5), close to Scotland and in Lake Baikal (see |301H303] for overviews). Another 
project, SAUND, has been using a very large but extremely sparsely instrumented hy- 
drophone array of the US Navy, close to the Bahamas |304H306] . The array of hydrophones 
covers an area of 250 km^, has good sensitivity at 1 - 500 kHz and can trigger on events above 
100 EeV with a tolerable background rate. 

An extended test configuration with various hydrophones, named AMADEUS, has been 
deployed together with the ANTARES detector |307| . The sensors are arranged in local 
clusters (size scale one metre), with inter-cluster distances between 12.5m and 340 m. First 
studies focussed on the ambient noise levels (see Fig 49 ) and their dependence on sea state 



and precipitation, with the conclusion that under calm conditions the noise in the relevant 
frequency range is below typical acoustic neutrino signals. In addition, the 3-dimensional 
arrangement of hydrophones allows for locating transient signals and investigating the rate 
of background signals of the expected bipolar shape that come from an appropriately defined 
fiducial volume. For the first time in the deep sea, the intrinsic background to acoustic neu- 
trino signals can thus be studied down to an amplitude level of a few 10 mbar, corresponding 
to detection thresholds at the EeV scale. 



N 



> 

(D 

m 

Q 

w 

Q. 



-70 
-80 
-90 
-100 
-110 
-120 
-130, 







' ' ' ! ' ' ' 


\ ' ' '■— 




































"bSS 


















Lv.,.i, 










< . 1 . < , 





I 



10^ 



^ ^10^ 



10" 



10 



20 40 60 

Frequency (kHz) 



80 



100 



Figure 49. Power spectral density 
(PSD) of the ambient noise 
recorded in AMADEUS. Shown in 
shades of grey is the occurrence 
rate in arbitrary units, where dark 
colours indicate higher rates. The 
white dotted line is the median 
value of the in-situ PSD and the 
black solid line indicates the noise 
level recorded in the laboratory 
prior to deployment. 
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8. Summary and Outlook 

Five decades after first proposals, and three decades after first practical attempts to 
build high-energy neutrino telescopes, we may be close to a turning point. IceCube, the first 
cubic-kilometre neutrino telescope, has started data taking in its final configuration, and 
preparations for counterparts on the Northern hemisphere are advanced. 

The strong case for high-energy neutrino astronomy has remained unchanged over all the 
time, but the requirements on the necessary sensitivity have tightened continuously. The 
detection of first high-energy neutrino sources from outer space lays still ahead and may, with 
some optimism, be expected in the next few years. Galactic "Pevatrons" as those detected by 
Milagro are in reach after a few years of IceCube data taking if the corresponding predictions 
are correct. Models assigning the most energetic cosmic rays to Gamma Ray Bursts will be 
tested within a couple of years. However, clear observations are all but guaranteed. 

Discoveries or non-discoveries of IceCube will have a strong impact on the future of the 
field and possibly mark a "moment of truth". Clear evidence for neutrino sources would 
pave the straight way for KM3NeT in the Mediterranean Sea, with possible modifications 
to the present design (high energies vs. low energies) according to the IceCube results. 
KM3NeT would then complement IceCube in its field of view; for results not depending on 
the hemisphere (e.g. diffuse fluxes, exotic particles) it could conflrm and extend the IceCube 
findings using a detector with different systematic uncertainties, and in addition it would 
explore the the central parts of the Galactic plane with superior sensitivity. 

Missing or marginal evidence for sources from IceCube may lead to different develop- 
ments. One option would be to continue the venue of detectors which explore the energy 
range characteristic for Gamma Ray Bursts and Active Galactic Nuclei. This would require 
configurations with significantly larger spacing than the present KM3NeT design, resulting 
in a penalty at energies below a few tens of TeV, and in sacrificing many Galactic source 
candidates. Abundance and characteristics of extragalactic sources are not expected to dif- 
fer much between the Northern and Southern skies. A factor of three or so in sensitivity 
compared to IceCube would therefore only provide a limited additional discovery window 
for extragalactic sources. Since firm fiux predictions are absent, the resulting "gambling" 
would be justified only if indeed a substantially larger sensitivity increase with respect to 
IceCube could be reached. 

Predictions for Galactic sources are much more firm, provided we assume dominantly 
hadronic emission. In contrast to extragalactic sources, we have good reasons to assume 
that we may be close to the discovery region, so that a sensitivity increase by a factor 3-5, 
with a telescope observing the Southern hemisphere, indeed counts. Therefore a second 
option would be to focus on Galactic sources like Supernova Remnants and tailor the array 
to these sources, i.e. keeping high sensitivity down to about a TeV. The uncertainty here 
comes from the question to what extent the observed gamma rays are indeed due to hadronic 
production. 

The third option wouki be an even larger leap in size than for the first option. It would 
address energies above 100 PcV with the help of new technologies like radio or acoustic 
detection and envisage 100-1000 cubic kilometres of instrumented volume. This option 
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might still have sensitivity to neutrinos from AGN jets but would also well cover the energy 
range of cosmogenic neutrinos. In contrast to optical detectors, new-technology detectors 
are still in the R&D stage and also have no natural calibration source like atmospheric 
neutrinos for optical detectors. 

The fourth option would define, at least for the time being, an end to the search for 
neutrinos from cosmic accelerators. The field would focus on optical neutrino detection with 
dense detectors optimised for investigating oscillations of accelerator neutrinos (Mediter- 
ranean Sea) and atmospheric neutrinos, or, even more pretentious, to study supernova bursts 
beyond our own Galaxy or even proton decay. 

For the moment, however, we don't see a reason to be pessimistic. We have made a 
factor-of-thousand step in sensitivity compared to a dozen years ago. This is far more than 
the traditional factor of ten which so often led to the discovery of new phenomena. For 
instance, looking back to our own field, the prospects for discovery had not been estimated 
overly high before launching the first X-ray rocket in 1962, or before detecting the Crab in 
TeV gamma rays in 1989. History told a different story, as we know by today. 

We have a good chance to open a new window to the Universe, but we don't know it for 
sure - and that, at the end, may be the most tempting and challenging situation that we 
possibly can imagine! 
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